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Abstract Ground-based 7-ray astronomy, which provides access to the TeV energy range, is 
a young and rapidly developing discipline. Recent discoveries in this waveband have important 
consequences for a wide range of topics in astrophysics and astroparticle physics. This article 
is an attempt to review the experimental status of this field and to provide the basic formulae 
and concepts required to begin the interpretation of TeV observations. 
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1 Introduction 

Ground-based 7-ray astronomy effectively began in 1989, with the first robust 
detection of a TeV 7-ray source, the Crab Nebula, usin g the 10 m diain eter air- 
Cherenkov telescope of the Whipple Observatory fWe ekes et al. 19891 ). Nine- 
teen years on, Cherenkov telescopes have been used to detect ~80 very high 
energy (VHE; E > 100 GeV) sources, firmly establishing a new astronomi- 
cal domain, and the first discoveries of sou rces using air-shower particle de- 
tectors have been made ( Abdo et al. 2007"al ). The key advantage of ground- 



based instrumentation over satellite-based GeV instruments such as EGRET 



(jHartman et al. 19991) and the rec ently launched GLAST (now Fermi) large area 



telescope (LAT) (jThompson 2004 ) is collection area. The typical effective collec- 
tion area of a single Cherenkov telescope is 10^ m^, almost five orders of magnitude 
larger than can realistically be achieved via direct detection in space. So far, the 
Imaging Atmospheric-Cherenkov Telescope (lACT) technique has proven to be 
the most powerful approach in this energy regime, with sensitivity of ~1% of the 
flux of the Crab Nebula, angular resolution for single 7-rays of around 5', energy 
resolution of ^15%, and the ability to locate sources with precision down to 10". 
The limitations of the lACT technique in comparison to air-shower particle de- 
tectors include a limited field of view (~ 5°) and a relatively poor duty cycle, 
with about 1000 hours of useful observations obtainable per year. 

The current status of VHE 7-ray astronomy is best summarized by Figure [H 
showing the distribution of different types of VHE 7-ray source on the sky. In 
addition to a collection of extragalactic sources at high latitudes — identified 
as active galactic nuclei, a number of clearly Galactic sources line the Galactic 
equator. Their low latitudes imply kpc distances and most are extended with 
respect to the ~5' resolution of the instruments, implying emission region sizes 
of 0(10 pc). These sources include shell-type supernova remnants, pulsar wind 
nebulae such as the Crab Nebula, binary systems, the Galactic Center and a 
number of unidentified sources. 

This review aims to summarize the basic emission mechanisms suspected to 
be active in these objects, to provide a cookbook of useful formula for calcula- 
tions in this field and to give an overview of the principle astrophysical results 
of TeV astronomy in the last few years. This article does not attempt to pro- 
vide a complete list of VHE 7-ray sources and their properties — given the 
rapid progress in the field this task is better served by online databases such 
as Te VCat — but rather aims to illustrate key properties and mechanisms using 
selected objects. The reader is also referred to recent reviews covering related sub- 
ject Sj_such_asjlie^evie^ astrophysics with ground-based detectors 



by lAharonian et al. (2008hl ) including a more detailed discussion of air-shower 



physics, the review of supernova remnants at high energy by Revnolds (20081 ) 



the discussion of whet her accelerati on in supernova remnants can account for 



galactic cosmic rays by|H illas (2 0051). the review of c osmic-ray propagation and 



interactions by Strong. Mo skal enko &: Ptuskin f2007l ) and the review of pulsar 
wind nebulae by Gaensler Sz Slane (20061 ). 



2 Generation and propagation of VHE 7-rays 



In this section some of the basic characteristics of acceleration and radiation pro- 
cesses relevant to VHE astronomy are summarized, providing simple expressions, 
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approximations, and references to detailed work. 



2.1 Particle acceleration and propagation 

Except for possible production by top-down processes such as the decay of heavy 
particles, VHE 7-rays are produced only in the interactions of accelerated charged 
particles - nuclei or electron^- with ambient matter or radiation fields. The 7-ray 
production rate reflects the product of the densities of cosmic ray (CR) particles 
and 'targets', tracking particles as they propagate away from the acceleration 
site due to diffusion in magnetic fields or convection. 7-ray sources will therefore 
be extended objects, with sizes and shapes governed by particle flow speeds, 
distribution of targets, and possibly by the finite lifetime of particles due to 
interactions or radiative cooling. 7-ray sources which appear point-like, given 
the ~5' resolution of current instruments, are either very distant (extragalactic) , 
or associated with compact targets such as dense molecular clouds or stellar 
radiation fields. 

Unless convective flows or bulk motion dominate, particle transport is governed 
by scattering off inhomogeneous ambient magnetic fields, resulting in diffusive 
propagation, (r^) = 2Dt. The diffusion coefficient D is determined by the average 
strength B of the magnetic field and its degree of turbulence 6B on length- 
scales comparable to the gyration radius Rg. The slowest possible diffusion in an 
isotropic medium is achieved in the Bohm limit rj = 1, where the mean free path 
of particles is given by the gyration radius, resulting in D w rjRgc/3 or, with 
Rg,pc ~ 0.0011i?Tev/-S_L_^G) ^nd particles of unit charge 

(r2jV2 ^ 0.01(r/i?TeVV/^MG)'/' (1) 

Typical interstellar fields ar e of the order of 3 p,G. The coeffic ient rj can be 
estimated as r/ ~ {dB/By^ (IStrong. Moskalenko fc Ptuskin 2nn7l ^. Bohm diffu- 
sion is very slow, compared for example to the speeds of young supernova shocks. 
Large-scale propagation in the Milky Way fe.g. lStrong. Moskalenko k, Ptuskin 2nn7l )| 



is governed by larger diffusion coefficients; in leaky box models one obtains 
D ~ A;102^£;ggv cmVs where k is a coefficient of order unity and a = 0.3 — 0.6, 
or (r'pc)^''^ ~ (-E^gytyr)^/^, with the energy dependence reflecting the scale distri- 
bution of magnetic turbulence. 

Diffusive shock acceleration is most likely the principle mechanism behind the 
production of high energy particles, with supernova remnants (SNRs) believed to 
be one of th e main sources o f CRs. For a recent overview of shock acceleration 
in SNRs see Revnolds A supersonic flow, for example due to ejecta from 



a supernova explosion or a pulsar wind, terminates in a shock, balancing the 
pressure of the ambient medium. Best viewed in the rest frame of the shock, 
material is streaming into the shock from upstream at velocity u, is compressed 
by a compression factor r and flows away from the shock with a speed reduced by 
the shock compression ratio. High-energy particles scatter off turbulent magnetic 
fields on both sides on the shock and may diffusively cross the shock many times. 
Each time they cross the shock and are isotropized by scattering in the medium 
on the other side, they gain an energy of order /S.E/ E ~ u/c. In each cycle, there 
is a finite probability for escape with the downstream flow, naturally creating 
a power-law spectrum of accelerated particles, N{E) ~ E~^ with F ~ (r -|- 



^Here and in the following, 'electrons' stands for 'electrons and positrons' 
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2)/(r — 1), with r = 4 and F w 2 for shocks with high Mach number. Neglecting 
radiative coohng and adiabatic losses, the distribution of accelerated particles is 
uniform in the region downstream of the shock, where particles are swept away, 
and extends over a scale D/u into the upstream region. The acceleration rate 
(!/£') (d-E/dt) ~ /D is governed by the rate of shock crossings, determined by 
the diffusion coefficient. A small diffusion coefficient will keep particles near the 
shock and ensure rapid return across the shock front. The maximum energy is 
governed either by the finite lifetime of the shock, by synchrotron losses in case of 
accelerated electrons, or by the gyro-radius exceeding the shock size. Assuming 
Bohm diffusion, maximum proton energies achievable in supernova shocks of age 
1000 t-i yr are of order 

^max ~ n^tg^^cTeV (2) 

where is the shock speed in units of 1000 km s~^ (e.g. Revnolds (20081 )): for 
electrons, synchrotron losses limit the peak energy to 

E^ax ~ lOO^zg^^G^'TeV (3) 



See Zirakashvili &: Aharonian 



In shock acceleration in SNRs, a significant fraction of kinetic energy of the 
flow can be transferred to high-energy particles. Once the energy density in 
part icles is comparable t o that in the shock, nonlinear effects start to play a role 
(see Caprioli et al. 20081 . and references therein). The overall compression ratio 



r is increased beyond 4, but particles scattered upstream of the shock decelerate 
the inflowing material and generate a precursor, reducing the compression ratio 
at the main shock. Particles with gyro-radii that are small compared to the size 
of the precursor experience only the reduced compression ratio r, resulting in a 
steeper spectral index, whereas for the highest-energy particles r > 4. Nonlinear 
shock acceleration hence produces concave spectra with F > 2 at low energy and 
r somewhat below 2 at high energy. 

Both acceleration speed and maximum energy increase for small diffusion coef- 
ficients, i.e. from high magnetic fields and maximal turbulence of the field on all 
scales, resulting in r/ ~ 1. In this context, turbulent field amplification by stream- 
ing CRs has received increasing attention. Both resonan t and non-resonant i nsta- 
bilities of the magne tic field are driven by CR currents ( Lucek k. Bell (2000l ). see 



Caprioli et al. (20081 ) for further references). In SNRs, magnetic fields can be es- 



timated fro m the cooling length-scale of electrons produced at the shock front (see 
for example Berezhko. Ksenofontov &: Volk 20031 . Vink &: Laming 20031 ) , with re- 



sults consistent with the predicted oc nv? dependence, where n is the ambient 
density and u the shock velocity. Magnetic field amplification boosts the max- 
imum energy of accelerated protons, but reduces the maximum electron energy 
due to increased losses. 

2.2 Electronic origin of 7-rays 

Electron s produce high-energy radi ation primarily via the inverse Compton (IC) 
process ( Blumenthal &: Gould 197ol ). by up-scattering ambient cosmic microwave 



background (CMB), infrared, optical or, in special cases. X-ray, photons. 

Energy loss rates: in the Thomson limit, b = 4lE(.Et /m'^c^ <C 1 or 6 ~ 
15Se TeV-E'reV ^ li with E'g and Et denoting electron and target photon en- 
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ergies, respectively, the energy loss rate is given by 

dE/dt = (4/3)aTC7'C/rad (4) 

A black-body target radiation field can be approximated by setting Et to 2.8kT. 
At higher electron energies, around 300 TeV, 10 TeV and 30 GeV for scattering 
off CMB, IR from dust and visible light, respectively, the Klein-Nishina (KN) 
regime begins and the energy loss rate is reduced, with a log E^ dependence in the 
deep KN regime. In general the cooling time, r = Ee/{dEe/dt), for IC scattering 
is given by 

V = EJidEjdt) « 3.1 • 10^f/-;^,v,^-3i?-Tev/KN (5) 
where the KN suppression factor /kn can be parametrized as 

/kN ~ (1 + 6)-'-' ~ (1 + 40^e,TevA:rev)"'-' (6) 



for b < 10^ ( Moderski et al. 20051 . and Erratum). The synchrotron cooling time 



is given by the very similar expression: 

Ty, = EJidEjdt) « 1.3 • 10^S;J^-i^v (7) 

Radiation spectra: scattering of mono-energetic electrons on a blackbody dis- 
tribution of target photons of temperature T results in a broad spectral energy 
distribution (SED) of the resulting 7-rays (Figure [2]). In the Thomson regime, 
energy losses are gradual and the 7-ray SED peaks at 

E^^TeY ~ SSE'g rp^v^^eV (8) 

resulting in 

Se,TeV ~ ll^'^TcV (9) 

for scattering of CMB photons. In the KN regime, electrons tend to lose a large 
fraction of their energy in a single IC event and the corresponding SED peak 
energy is shifted to: 

E^,TeV ~ ^e,TeV^^^^2.^16y).8)l/0.8 ^^^^ 
The equivalent expression for the typical energy of synchrotron photons is 

^s,eV = 0M7Elr^,yB^G (11) 

Whilst the energy loss rate of electrons depends (in the Thompson regime) only on 
the energy density in target radiation fields, the spectrum of 7-rays is strongly 
influenced by the spectral distribution of target photons (Figure E]). Since in 
the Thomson limit the 7-ray energy scales with the square of the electron en- 
ergy, the IC 7-ray spectrum is harder than the spectrum of parent electrons: 
assuming isotropic distributions of electrons and target photons, a power-law 
distribution of electrons Ne{E) ~ E~^^ will generate IC (and synchrotron) spec- 

tra of index F^ = (Fg + l)/2, cc E-y . Deep in the KN regime, 7- 

ray spectra steepen by AF ,^ « (Fe + l)/2, ~ E.y^^''~^^\ up to logE'^ terms 
(|Bhimenthal fc Gould 19701 '). Therefore, even for a pure power-law distribution 
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of electrons, 7-ray spectra exhibit a break corresponding to the transition to the 
KN regime. 

IC 7-rays of energy (in the Thomson regime) and synchrotron photons 
of energy -Egync probe the identical electron population provided that -E-y^TeV = 
380£'sync,ev(^^ev/-B^G)- For example, 11 TeV electrons produce both 1 TeV IC 
photons (on the CMB) and 1 keV synchrotron photons, if B=WO fiG. In this 
case, the ratio of energy flux in the synchrotron and IC bands is just C/mag/fA-ad) 
allowing the magnetic field strength in the source region to be determined. For 
the smaller (several //G) fields found typically in the ISM, X-ray synchrotron 
emission traces electrons of considerably higher energies than the TeV IC 7-rays 
and a model for the shape of the electron spectrum is required to deduce the 
magnetic field. 

Which target photon fields are most relevant for 7-ray production depends on 
the electron spectrum and on the 7-ray energy considered. The relative yields 
of CMB target photons, IR target photons and visible target photons vary de- 
pending on the location of sources, e.g. with Galacto-centric radius, and with 
the proximity of strong radiation sources. At the location of the sun, typical 
radiation energy densities are 0.3eVcm~^ for IR from dust, and 0.3eVcm~^ for 
visible light, see IPorter. Moskalenko k Strong (2006l l (c.f. 0.26 eV cm for the 



CMB). For electron spectra without a cutoff, CMB photons will often dominate 
at sufficiently high energies, as the contributions of IR and visible photons are 
suppressed in the KN regime. In this case, 7-rays directly map the distribution 
of VHE electrons. For electron spectra with a cutoff, the situation can, how- 
ever, be reversed, with the efficient transfer of energy from an electron to a 7-ray 
over-compensating for the drop in scattering cross-section. 

Electrons passing through a medium containing atoms or plasma will also cre- 
ate 7-rays by bremsstrahlung ( Blumenthal &: Gould 1970l ). The loss timescale is 



energy-independent, with loss rates and spectra depending on the state of the 
medium (mainly the shielding length for ions). For neutral hydrogen atoms of 
density n per cm^, r = Ee/{dE(,/dt) ~ 3.9 • 10'^n~^ years. IC scattering on 
the CMB dominates in the 1 TeV emission of electrons in (neutral) media as 
long as n < 240 cm^'^. Bremsstrahlung 7-ray spectra produced by a power-law 
distribution of electrons are also power-laws, with F^ = Fg. 

2.3 Modelling radiation spectra from electron populations 

Energy losses by synchrotron radiation or IC will modify electron spectra com- 
pared to the index at injection, Fi, above energies where the energy loss time (Eq. 
[7]) becomes comparable to the age T of the electron source. If synchrotron losses 
or IC losses in the Thomson regime dominate, {dE /dt) sync+ic = —kE"^, the elec- 
tron energy distribution will be cut off (for burst injection in a source of age Tg) 
at E = 1/(kTs), the energy where the radiative lifetime {1/ E){dE / dt) equals the 
age Tg , or will exhibit a spectral break at the same energ y, with index T p increas- 



ing by 1 (in case of continuous injection over time Tg) (iKardashev 19 62). This 
will cause a cutoff or break with index change AF^ = 1/2 in IC and synchrotron 
spectra at the corresponding energies (equations [51 and 1 1 1 1) . The situation is more 
complicated if IC losses dominate over synchr otron losses, as might be the case 
near strong radiation sources, as discussed by Moderski et al. (20051 ). Entering 



the KN regime, electron energy losses are reduced, resulting in a harder "cooled" 
spectrum Fe f« Fj 1 AF with AF -1... - 1.5, and F^ « Fi. With IC 
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losses scaling slower than E'^, synchrotron losses will, however, always dominate 
above a certain energy, causing to change to f« Fi + 2. Figure [3] shows the ef- 
fects of both IC and synchrotron dominated cooling on the SED from an injected 
power-law of relativistic electrons. 

The above discussion assumes particle injection with a power-law spectrum, 
followed by radiative losses. If acceleration timescales and loss timescales are 
comparable, acceleration and losses cannot b e factored and complex spectra l 
patterns can arise, as discussed, for example, by Zirakashvili &: Aharonian (200?! ). 



In modelling radiation spectra of astrophysical sources, several approaches are 
followed: 

• One can determine the spectrum (and composition) of particles which, in- 
teracting with suitable target fields and radiation fiel ds, gives rise to the 
observed wide-band SED (e.g. lAharonian et al. 2005dl ). 



One can determine the energy spectrum of injected particles which, after 

accounting for the modification of the spectrum over tim e due to radiation 

losses and interactions etc., creates the observed SED. (e.g. Hinton &: Aharonian 2007l .| 
and in Figure [3|) . 

One can model the full dynamics of the sourc e, including particle acc elera- 
tion mechanisms, losses and interactions (e.g. Berezhko fc Volk 20061 ). 



The first approach is self-consistent but has the potential problem that, due 
to the many degrees of freedom in the choice of spectral parametrisation and 
targets, it may be difficult to arrive at a unique solution. Frequently, spectral 
parametrisation assumes a broken power law, with the fit sometimes resultin g in a 
large increase in spectral index at the break (see e.g. Aharonian et al. 2005dl ). It is 



non-trivial to see which mechanisms would produce such spectra; at synchrotron 
cooling breaks in electron spectra, the index increases only by one. 

The second intermediate approach attempts to cure this deficit by including 
(some of) the mechanisms which cause spectral breaks and cutoffs. However, 
simplifying assumptions are often made concerning the history of particle in- 
jection. In sources where particles are confined within an expanding envelope, 
such as SNRs or pulsar wind nebulae, adiabatic losses usually need to be taken 
into account, but frequently are not. Cooling during acceleration may mod- 
ify injection spectra compared to the often assumed power law. This approach 
usually gives reliable answers only if cooling times are short compared to the 
evolution/expansion timescales of the system, but long compared to acceleration 
timescales. 

The third approach is "best" in that it attempts to account for all relevant 
effects, but is most demanding and has the disadvantage that solutions are usually 
numerical and that it is non-trivial to understand the systems in terms of a clear 
one-to-one connection between input assumptions and predicted radiation SED. 

2.4 Hadronic origin of 7-rays 

An alternative source of VHE 7-rays are interactions of high-energy protons 
and nuclei with interstellar material. Through this mechanism, 7-ray emis- 
sion traces CR acceleration sites and CR propagation, as demonstrated in the 
GeV energy range by the study of diffuse 7-ray emission from the Milky Way 



(jHunter et al. 19971 ) . The interaction cross-section of VHE protons on hydrogen 
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nuclei of density n per cm^ is only weakly energy-dependent, app^inei ~35 mb, 
resulting in a lifetime of r w 3 • 10^n~^ yr at multi-TeV energies. The number 
of secondary particles produced per interaction increases with energy. Typically, 
half of the energy Ep of the primary is carried away by a leading nucleon, the other 
half is split between charged and neutral pions and a small fraction of heavier 
hadrons. This implies that about 1/6 of the primary energy is carried by a num- 
ber of 7-rays produced in vr^ decays. To a good approximation, and for E^ ^> m,r) 
the distribution in energy of 7-rays produced per interaction is scale invariant, 
dnry/dE^{E^, Ep) = D^{x = E^/Ep ). Parameterisations for the fragmen tation 
function D are given for example bv lKelner. Aharonian &: Bugavov (20061 ). Fig- 
ure [2] illustrates the SED of 7-rays produced by mono-energetic protons. The 
SED of secondary 7-rays peaks at about 1/10 of the energy of the primary. For a 
power-law proton distribution dup/dEp = kpE~^ the 7-ray distribution is again 
a power law with the same index: 



= cnapp^inei I ( ^ ] = cnapp^ineikpE^ [ D{x)x^ 

E^at JE-y CLEp yEp J Ep 'Jo 

(12) 



It is convenient to write this as 



dn^ , / E, 



CnCTpp^inell^kp ( j (13) 



dE^dt y f 

since for / ~ 0.15 the coefficient k ~ 6 is approximately constant for F between 2 
and 4. This property should not be misinterpreted in the sense of a delta-function 
approximation, namely that in general the 7-ra y spectrum at energy E^ traces the 
proton spectrum at energy Ep/ f. In particular, Kelner. Aharonian &: Bugavov (2006l )| 



show that an exponential cut-off in the proton spectrum, exp{—Ep/Ecut) is trans- 
formed into a more gradual cut-off, exp{—{16E^ / Ecut)^^^) in the 7-ray spectrum. 
Obviously, given the width of the 7-ray SED illustrated in Figure El any feature 
in the proton spectrum will re-appear smoothed in the 7-ray spectrum. 

One application is the estimation of the 7-ray flux from molecular clouds illu- 
minated by CRs, possibly from a nearby accelerator. For a cloud of mass M and 
distance d illuminated by a proton flux ^{E) = kpE~^ (per area, solid angle and 
energy), one obtains a 7-ray flux 



^,{E) = ^-^^Kkp ^ (14) 



with the hydrogen mass rriH- With a small correction for heavier nuclei and 
assuming a CR proton flux as measured on Earth, this translates, for example, 
into a flux 

^>^(> ITeV) ^ 1.6 • lO^^^Me/dkpc cm'^s"^ (15) 
where Mg is the cloud mass in units of 10^ solar masses (see e.g. Aharonian 199ll ). 



These estimates assume that the CRs permeate the cloud; the increased magnetic 
fields inside dense clouds may however influence the diffusion coefficient and hence 
the spectrum and flux of CRs. In the case of clouds illuminated by nearby 
CR accelerators which were active for a relatively short time, energy-dependent 
diffusion will modify the spectrum at the cloud as compared to the spectrum 
generated at the accelerator; high-energy particles reach the cloud first, resulting 
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in a harder 7-ray spectrum at the cloud (see iGabici. Aharonian fc Blasi 20071 . 
and references therein). 

Dense environments wih also result in production of bremsstrahlung 7-rays by 
electrons; assuming that electron and proton spectra are identical and that the 
ratio of electrons to protons is kep, one obtains a flux ratio $7,brems/^7,hadr ~ 
3.3 kep at 1 TeV, with only weak dependence on photon energy due to the changing 
CTpp (in the region far from cut-offs). 

2.5 Absorption during propagation and opacity of sources 

The mean free path of VHE 7-rays in hydrogen gas is governed by the electron- 
positron pair production cross-section and has a value of 80gcm~^ or equiv- 
alently 5 • 10^^ hydrogen atoms cm~^; for all practical purposes, the universe 
is transparent to this process. The more relevant process is the absorption 
by pair production on ambient (CMB, IR, visible or X-ray) photons of en- 
ergy Et- The process acts above the threshold E^Et = m?c^ or E^^TeV = 
0.26/ Er. f^v and the absorption cross section for an isotropic photon field a-y^ 



(see e.g. lAharonian. Khangulvan &: Costamante 20081 ) peaks close to threshold 



at Ej^TeV = 0.9 /ET,eV- Figure m illustrates absorption by different (blackbody 
and power-law) photon spectra. 

Two situations where absorption is important are (a) the observation of extra- 
galactic sources and (b) TeV emission near intense sources of radiation. Figured] 
also shows the optical depth for VHE 7-rays interacting with intergalactic radi- 
ation fields taking into account cosmological evolution of the background fields; 
the range of 7-rays is about z = 0.03, ~ 0.1, ~ 1 for E^ = 10, ~ 1, ~ 0.1 TeV. At 
PeV energies, the range is reduced to galactic distance scales. At energies up to 
tens of TeV, on the other hand, absorption of galactic sources is almost negligible 
(Figur e HI), even for sources near the Gala. ctic Center with its increased radiation 



fields (jMoskalenko. Porter fc Strong; 200d ). 



Absorption is relevant for 7-rays produced in radiation-intensive environments, 
for example in systems of compact objects in close orbit around massive stars. 
Radiation produced, for example, in a massive X-ray binary at 0.1 AU from 
the massive star traverses 0(10^^) eVcm~^, resulting in a large optical depth 
at TeV energies if the 7-ray source is behind the star, allowing head-on colli- 
sions between 7-rays and stellar photons. Because of the peaked absorption cross 
section, narrow (black-body) background photon spectra cause selective absorp- 
tion in part of th e VHE ener g y range and lead to significantly modified spectra 
(see Figure HlandlAh aronian. Khangulvan &: Costamante (20081 )). In such dense 
absorbers, the electrons produced may undergo IC scattering, creating a pair 
cascade which proceeds until 7-ray energies are low enough that the medium 
becomes transparent. The absorption dip in the 7-ray spectrum is then ac- 
companied by a correspon ding enhancement at lower energies (see for example 
Protheroe &: Stanev 19931 ). Absorption within the source is also a key factor for 



7-ray emission from compact regions (knots, blobs...) within AGN jets, due to 
the high radiation density in these regions; see Section] 



3 Instruments for TeV astronomy and their characteristics 

Given the very low fluxes of 7-rays in the VHE regime — 0(10"") photons per 
cm^-second (a few photons per m^-year) above 1 TeV for strong sources, direct de- 
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tection by space-based instruments is excluded. Ground-based instruments detect 
secondary products resulting from the development of 7-ray initiated air-showers; 
either particles reaching the ground or Cherenkov light emitted by shower par- 
ticles in the atmosphere. In contrast to the well-collimated electromagnetic air- 
showers induced by 7-rays (or electrons), air-showers initiated by CR nucleons 
typically feature a number of electromagnetic sub-showers induced by vr^ decays 
and contain muons from charged pion decays (see Figure [6]). Rejection of the 
background of showers initiated by charged CRs is a key performance criterion 
for 7-ray detection systems, and is usually achieved on the basis of shower shape 
or muon content. A more detailed discussion of air-s hower characteristics and the 
detection systems used can be found for example in Aharonian et al. (2008hl ). 



3.1 Instrument characteristics 

For ground-based instruments detecting 7-rays via their shower development in 
the atmosphere, effective detection areas, A(E) (defined such that the differential 
detection rate R^{E) = ^ry[E)A{E)), have a sub-threshold region where they 
exhibit a steep rise with energy, and a high-energy region where A{E) varies 
only weakly with energy. In the sub-threshold region, the detector triggers only 
because of favorable fluctuations in the development of an air-shower. In the 
high-energy region, every air-shower within a certain fiducial region is recorded. 
The 'energy threshold' of a detection system is usually quoted as the energy at 
which the peak detection rate R{E) occurs for typical power-law 7-ray spectra. 
The threshold thus determined obviously depends on the assumed spectral index, 
but is always located in the transition region between the steeply rising part and 
the nearly constant region of A{E). Individual events may be detected at energies 
well below this nominal threshold. 

Two criteria govern detectability of a source during an exposure time T: (a) a 
minimum number no (usually 5. ..10) of 7-rays must be detected, TR^ > no, and 
(b) the 7-ray signal must be be significant above fluctuations in the background, 
approximately (TR^) / (TrjcRRcR^)^^'^ > ^o- Here, -RcR is the detection rate of 
background CRs per solid angle, tjcr the efficiency for CRs passing analysis cuts 
relative to the corresponding efficiency for 7-rays, and 0,, the solid angle over 
which the signal from a source has to be integrated, given either by the point 
spread function (PSF) of the instrument or the source size. For point sources 
and a Gaussian PSF, O ~ ^r^is where 6qs is the 68% containment radius of the 
PSF. Current instruments typically operate in the background dominated regime, 
implying that minimal detectable fluxes scale as T~^^'^riQ^6Qs- For sources which 
are large compared to the PSF, J7 ~ vr^^, and sensitivity degrades linearly with 
source radius Og. 



3.2 Cherenkov imaging of air-showers 

In the past decade. Imaging Atmospheric Cherenkov Telescopes (lACTs) have 
emerged as the most powerful instrument for pointed observations of VHE 7- 
ray sources. lACTs use focusing mirrors to image the Cherenkov light emitted 
by shower particles onto a pixelated photon detection system (see Figure E]) ; a 
summary of the characteristics of current and selected previous instruments is 
given in Table [TJ 

At the maximum of the shower development, around 10 km a.s.l. for TeV ener- 
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gies, the Cherenkov threshold for electrons is around 40MeV and the Cherenkov 
angle is 0.7° or less. Light emitted at the Cherenkov angle reaches the ground 
within a circle of 100 to 150 m depending on the height above sea level of the de- 
tection system. Multiple-scattering angles of shower particles near the Cherenkov 
threshold are comparable to the Cherenkov angle, resulting in a more or less uni- 
formly filled light pool, with typically 10 detected Cherenkov photons per TeV 
shower energy and m^ mirror area for photomultiplier sensors. With increasing 
energy, the central density in the light pool is enhanced due to deeper penetration 
of showers. Triggering and image reconstruction usually requires 50 to 100 de- 
tected photons and sets the scale for the dish size. The pixel size of the detection 
system should be matched to the size of features in air-shower images; simulation 
studies show saturation of performance for pixels much below 0.1° diameter — 
close to the typical rms width of a 7-ray image at TeV energies. The asymptotic 
collection area for lACTs is determined by the maximum impact distance for 
which shower images still fall within a camera and hence by the camera field of 
view (FOV). At 2000 m a.s.l. the impact distance limitation is approximatey 
100 m per degree of the opening angle of the camera field of view (for showers 
close to zenith). 

Most modern instruments use multiple telescopes (a) to image the air-shower 
from different viewing angles for improved reconstruction of 7-ray direction and 
rejection of CR background and (b) to apply a coincidence requirement rejecting 
single-telescope triggers caused by CR muons with impact points close to a tele- 
scope mirror, or by night sky background. Telescope spacing needs to be large 
enough to provide a sufficient baseline for stereoscopic measurements, but small 
enough that multiple telescopes fit within the Cherenkov light pool; the exact 
spacing tends to be uncritical within a range of ~ 70 m to 150 m. Depending on 
selection cuts, telescope systems such as H.E.S.S. (Tabled]) provide an angular 
resolution for single 7-rays of 3' to 6', a 7-ray energy resolution of around 15% and 
a CR rejection factor of 0(10"^) or better. Combined with the energy threshold 
of about 100 GeV and a high-energy effective area of some 10^ m^, this allows de- 
tection of sources of 1% of the strength of the Crab Nebula {uF^ ~ 3 x 10^^^ erg 
cm~^ s~^ around 1 TeV) within 25 h of observations close to Zenith. Performance 
of lACTs is at some point limited by fluctuations in air-shower development; for 
example, at energies below 10 GeV an air-shower is in principle still detectable, 
but because of the small number of Cherenkov-emitting particles, the energy de- 
termination is quite unreliable. Current instruments, however, are still relatively 
far from reaching fundamental limits; for example, shower fluctuations allow an 
angular resolution of a fraction of an arc-minute at TeV energies, provided that 
the number of Cherenkov photons d etected is suffic ient, and a CR rejection fac- 
tor of order 10~^ at TeV energies (jHofmann 20061 ). A non-trivial issue in the 
analysis of data is the absolute energy calibration. Cherenkov light from local 
muons is often used to calibrate the response of instruments, but variations in 
atmospheric profile and transmission and of the orientation of the shower axis 
relative to the geomagnetic fie ld can influence the shower development and the 
light yield fe.g. iBernlohr 2000l ) and cause systematic calibration uncertainties at 
the 10-20% level. 
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3.3 Ground-based particle detectors 

The direct detection of air-shower particles offers a method of 7-ray detection with 
close to 100% duty cycle and very wide (~1 steradian) field of view. Because of 
these advantages this method offers an interesting complementary approach, de- 
spite the fact that the point-source sensitivity of such detectors is currently almost 
two orders of magnitude poorer than the most sensitive lACTs (at a few TeV). 
Indeed the extensive air-shower sampling (EAS) techniq ue has recently prod uced 
its first contributions to the catalogue of TeV sources (|Abdo et al. 2007al ^ and 



proved its usefulnes s for all sky sur veys (at modest sensitivity) and the detection 



of diffuse emission (lAbdo et al. 2008 ) 



The main challenges of the EAS approach are the rejection of the CR back- 
ground and directional and energy reconstruction using the exponentially decreas- 
ing tail of particles detected well beyond shower maximum. High altitudes are 
therefore critical to achieve low (< 1 TeV) thresholds with such instruments. In 
addition, whilst in principle all-sky detectors, the field of view obtained is limited 
in practice by the rapid increase in energy threshold wi th zenith angle; typ ically 
a factor two increase between 0° and 30° zenith (see e.g lAbdo et al. 2007al ). The 



reconstruction of the primary 7-ray direction is based on shower front timing. 
The arrival time of the shower front can be determined with an accuracy of a few 
nanosec onds over ~ 100 m baselines leading to typical resolutions of 0.3° - 1° 
(see e.g. I Atkins et al. 20031 ^. The much larger fluctuations present in the particle 



number at ground level with respect to the essentially calorimetric air-Cherenkov 
approach, make primary energy determination extremely difficult for EAS detec- 
tors. The rejection of the hadronic background is based on the muon content of 
showers and/or the distribution of shower particles on the ground. 

The Water Cherenkov approach pioneered by the MILAGRO collaboration, ap- 
pears to represent the most cost effective method of achieving complete ground 
coverage. MILAGRO is a 80 x 60 x 8 m pond instrumented with PMTs and sur - 
rounded by 175 water tanks, located at an altitude of 2630 m ( Atkins et al. 20031 ). 



MILAGRO achieves its best background rejection power and sensitivity in the 
regime above 10 TeV. Widely spaced (~1 m^ detectors > 5 m apart) scintillation- 
based detectors have been used for CR measurements and in the search for UHE 
7-ray sources for many years Instruments of this type located at high altitudes 
can be used for 7-ray as tronomy around a few TeV, as demonstrated with the 
Tibet Air-Shower Array ( Amenomori et al. 19991 ) at 4300 m above sea level. The 



ARGO-YBJ instrument is a 5800 m^ complete ground coverage instrument at the 
same site, utilizing resistive plate counters (RPCs) and achieving a threshold of 
a few hundred GeV and a sensitivity similar to that of MILAGRO. However, this 
approach is likely prohibitively expensive for a much larger area next generation 
detector. 



4 VHE 7-ray sources 

As can be seen from Figure [T^, the TeV sky, despite a modest number of known 
objects (~80), contains a diverse collection of different object classes. Numeri- 
cally dominant are the sources clustered along the plane of our galaxy (see Fig- 
ure [TJ)). In contrast to GeV energies, sources, rather than diffuse emission, domi- 
nate our current view of the Galaxy. In the following, we will attempt to summa- 
rize the important characteristics of each class of TeV 7-ray sources in the context 
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of the discussion on particle acceleration, transport and radiation given in Section 
El Table [2] gives a summary of prominent Galactic VHE 7-ray emitters for which 
firm identifications exist. There is insufficient space here for an adequate dis cus- 
sion of diffuse TeV emission. We refer the reader to ^ Aharonian et al~ for 
a review of this topic and also the complex situation at the Galactic center. 



4.1 Supernova remnants 

Ever since IZwickv (19^ ). supernova remnants have been viewed as the most 
likely sources of galactic cosmic rays up to an energy of at least that of the knee 
of the CR spectrum around 10^^ eV, and possibly beyond 10^ eV. This argu- 
ment is based in part on the energy input required to maintain the CR flux in 
the Galaxy, dE/dt ^ pVr ~ 5 • 10^°ergs"\ with the energy density in cosmic 
rays p ~ 1 eV cm~^, the confinement volume of cosmic rays V ~ 10^^ cm^ and the 
characteristic residence time r ~ 10'' yr. With a kinetic energy output of 10^^ erg 
per explosion and a rate of a 2-3 per 100 years, an average of 10% of the super- 
nova kinetic energy needs to be converted into CR energy. The first-order Fermi 
acceleration process outlined in Section 12.11 possibly enhanced by magnetic field 
amplification, provides a means to plausibly reach this efficiency and a maximum 
energy beyond 10^^ eV, and naturally provides a power-law spectrum with an in- 
dex around 2, which can explain the observed CR index if energy-depende nt diffu- 
sion and escape is invoked (see e.g. IStrong. Moskalenko fc Ptuskin 200?! ). While 
this scheme is generally accepted, there are still open questions concerning the 
consistency of required diffusion speeds and (the lack of a significant) anisotropy 
of the CRs at the Earth. Equally important, while the presence of high-energy 
electrons in SNRs can be inferred from the non-thermal X-ray spectra measured 
in several objects, modulo the a priori unknown, strength of -B-fields, efficient 
acceleration of nuclei still lacks undisputed evidence, and observational proof is 
missing that acceleration in SNR can quantitatively account for the observed CR 
spectrum. 

VHE 7-rays trace the relevant populations of energetic particles in SNRs. The 
flux of hadronic 7-rays is proportional to the CR density times the density of tar- 
get gas, and the electronic IC 7-rays directly trace the electron density, given that 
target photon fields are likely to be almost constant on the scale of the SNR. The 
ratio of IC 7-rays and synchrotron X-rays is determined by the strength of the 
magnetic field. A detailed study of acceleration in an SNR shell, and distinction 
between the shell and the nebula of a possible pulsar created in the explosion, re- 
quires the shell to be resolved in VHE 7-rays. This has s o far been achieved 
for four SNRs: R X J1713.7-3946 (lAharonian et al. 2006bl^. RX J0852.0-4622 
(alias Vela Junior) (jAharonian et al. 2007dl'). RCW86 (lAharonian fc et al. 20081 ) 
and, most recently, SN 1006 ( Naumann-Godo &: et al. 20061 ). The 7-ray spec- 
trum measured for RX J1713.7— 3946 extends to several tens of TeV and follows 
a power law with index 2 up to about 20 TeV where a cut-off sets in (Figure 
[71), demonstrating the p resence of primaries in the 100s of TeV energy range 
( Aharonian et al. 2007el ). In all cases, there is a strong similarity between the 
morphology observed in VHE 7-rays and in non-thermal X-rays, once the dif- 
ference in PSF of the instruments is taken into account (see Figure E]). The 
correlation between X-rays and 7-rays seems to argue in favor of a leptonic origin 
of 7-rays. Correlation with gas density, as traced by CO, is less pronounced. 
This may, however, be due to a lack of resolution along the line of sight (i.e. in 
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distance estimated by line velocity), making it difficult to know if a given gas 
mass is actually co-located with the accelerated particles. Only in a few special 
cases can the association between gas and remnant be proven by an observed 
high velocity dispersion of the gas, resulting from the interaction with the SNR 
shock, or by OH maser emission. 

The complex morphology of the remnants, with non-uniform emission along the 
rim, reflects their interaction with their environment and makes interpretation 
difficult, in particular given the lack of high-resolution, 3-dimensional information 
on the surrounding gas density. In the case of SN 1006, located off the Galactic 
plane in a less complex environment, the observed bipolar emission pattern can be 
modeled more easily: SN 1006 is believed to be threaded by a relatively uniform 
magnetic field, and non-thermal emission marks the polar caps where the B-field 
vector is roughly parallel to the expansion direction. In the equatorial regions, 
where field lines are perpendicular to the direction of shock propagation, particle 
injection into the Fermi process is presumably inefficient, since particles spiral 
along field lines and are immediately overrun by the shock. In the polar regions, 
injection is more effective, resulting in significant CR current and most likely in 
turbulent magnetic field amplification, rendering the acceleration process even 
more efficient. 

Explaining the VHE emission from RX J1713.7— 3946 — probably the best 
studied case — as IC emission of electrons requires magnetic fields of around 
10 /xG ( Aharonian et al. 2006bl ). These magnetic fields are at variance with the 

significantly higher fields of at least 50 fiG determined from th e width of X-ray 

emitting filaments, translated into an electron cooling time (e.g. Berezhko &: Volk 20061 ) .1 
and from the temporal variations observed in X-ray structures ( Uchivama et al. 2007l ).| 
which directly measure the cooling time. Questions remain, though, whether 
these high magnetic fields are characteristic of the entire SNR, and if the ob- 
served structures are due to variations in electron density rather than magnetic 
fields; only if this is the case, an IC origin of the 7-rays is firmly ruled out. There 
are no significant variations of 7-ray spectral shape across the TeV SNRs and the 
observed spectra are well described by hadron acceleration models which gener- 
ate a power-law distribution with an index close to 2, with a gradual high-energy 
cut-off (Figure [7]) . Leptonic models tend to generate spectra that are too hard at 
low energy, reflecting the fact that the index of IC 7-rays is 1.5, for an electron 
index of 2 (see Section [2.2p . This problem can be solved to a limited extent by 
adjusting the composition of target radiation fields and by introducing multi- 
zone models with different cut- offs, where different IC peaks are su perimposed to 
mimic an E~'^ spectrum (e.g. Porter. Moskalenko &: Strong 20061 ). In hadronic 
models, on the other hand, the strong X-ray /7-ray correlation is non-trivial to 
obtain. In the case of high B-fields, electron lifetimes are comparable to accel- 
eration timescales and the X-ray flux is influenced both by the energy input in 
accelerated electrons and by the strength of B-fields. Hadronic 7-rays, on the 
other hand, refiect proton fiux — which should scale with the (injected) electron 
flux — multiplied by the gas density. A strong correlation between X-rays and 
7-rays requires a link between gas density and B-field strength; magnetic field 
amplification (see Section [2.ip may provide a mechanism for this. 

A final demonstration of CR origin in SNRs may be achieved by a combina- 
tion of wider spectral 7-ray coverage, improved resolution or morphology, and 
larger-scale measurements of magnetic fields. Detection of neutrinos from SNRs 
would also demonstrate hadronic origin, but is challenging even for the largest 
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instruments such as ICECUBE, and suffers similarly from the fact that, for a 
quantitative analysis, the target gas density needs to be known. Investigation 
of global CR energetics and spectra will in any case certainly require the 7-ray 
detection and spectral analysis of a representative sample of SNRs. 

Another approach towards demonstrating CR acceleration in SNRs is to look 
for dense molecular clouds adjacent to, or interacting with, an SNR. In clouds, 
interactions of accelerated protons and nuclei will give rise to an enhanced 7-ray 
flux proportional to the cloud's mass (Equation [T5|) whereas IC radiation from 
electrons i s not enhanced. Two ca ndidate syster ns where this m ight be occurring 
are W28 (|Aharonian et al. 2008el l and IC443 (|Albert et al. 2007ai ) (Figure ED . 



W28 is an old remnant (30-150 kyr) which has most likely released most of its 
CRs. VHE 7-ray data show four emission hot-spots coincident with enhancements 
of gas density; if interpreted as proton interactions in passive clouds, their masses 
imply a CR flux which is 10 to 30 times the flux near Earth, a plausible value 
given the proximity (at least in projection) of the remnant. A similar situation 
is seen in IC 443, where TeV emission coincides with a massive molecular cloud, 
with OH maser emission indicating that the SNR shock wave is hitting the cloud. 

4.2 Pulsars and pulsar wind nebulae 

The first VHE 7-ray source to be detected, the Crab Nebula ( Weekes et al. 19891 ). 



is a pulsar wind nebula (PWN) , where populations of elect rons and positions with 
energies up to PeV energies emit X-rays and IC 7-rays (e.g. Atovan &: Aharonian 19961 ) .[ 



As PWNe have a well-defined central energy source and are typically close enough 

to be spatially resolved, they allow relativistic fiows, and the s hocks which result 

when these winds collide with their surroundings, to be studied ( Gaensler Slane 20061 ) .[ 



PWNe are the most abundant class amongst the sources discovered in the H.E.S.S.| 

survey of the Galactic Plane. Many of the fundamental concept s co ncerning 

PWNe are summarized in t he seminal papers bvlRees fc Gunn ^19741 ) and lKennel &: Coroniti ri984l ) 



for a recent summaries see iGaensler fc Slane (2m(i )'. 



A supernova explosion may create a pulsar, a neutron star with a magnetic field 
axis which is misaligned with the rotation axis. The rotating magnetic dipole will 
emit electromagnetic radiation at a luminosity S ~ 3 • W"^^ Bi2P^t ^rg and 
will spin down; B12 is the surface magnetic field in units of 10^^ G and Pms the 
period in ms. Pulsars with 7-ray PWN tend to have E around and above 10^^ 
gj.gg-i_ "Yhe time dependence of the spin-down energy loss is given by 

E{t) = Eo/{l + t/Tr (16) 

with p = {n + — 1), with the characteristic spin-down time r. For pure 

dipole radiation one has a breaking index n = 2 and r ~ 10Pq^^/Bi2 y where 
-fb,ms is the birth period of the pulsar; measured values for n lie between 2 and 
3. The rotating field creates a voltage drop of order 10^"^ B12 / P^^. V which can be 
used to accelerate particles, fed by electron-positron pair cascades in the giant 
electric and magnetic fields near the pulsar surface. Accelerated electrons will 
emit 7-radiation which appears pulsed to a stationary observer away from the 
rotation axis. Both the polar caps of pulsars and the "outer ga p" at a couple 
of stellar radii have been considered as emission regions (see e.g. Harding 200?! . 



for a review and references). Radiation emitted near the polar cap of the pulsar 
should cut off sharply at a few GeV since cascading in the high fields in this 
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region prevents escape of higher-energy photons. In the region of the 'outer gap', 
energies of tens of GeV may be reached. 

The relativistic electron-positron wind from the pulsar terminates in a shock 
where the ram pressure of the wind is balanced by the pressure of the surrounding 
nebula. At the shock, the kinetic energy of the wind is transformed into random 
motion. Outside the shock, the resulting relativistic electron-positron gas will 
convect outwards at subsonic speeds, v < c/\/3, v decreasing initially as 1/r^, 
forming an expanding PWN visible in sychrotron radiation and IC 7-rays. As- 
suming Bohm diffusion, convection of pa rticles will dominate over diffusiv e prop- 
agation in most of the nebula (see e.g. de Jager &: Djannati-Atai' (20081 )). The 
reverse shock created in the expanding SNR may collide with the expanding PWN 
after some kyr and may temporarily halt t he expansion of the PWN. The evol u- 
tion of P WNe is summarized for example in lBlondin. Chevalier &: Frierson (200 1 ) .| 

In retrospect, it appears plausible that PWN are dominant among galactic 
VHE 7-ray sources: the energy content of a PWN - of order 10^^ erg - is small 
compared to the energy of ~ 10^^ erg dissipated in a supernova shock, but since 
a large fraction of the energy is carried by relativistic electrons with radiative 
lifetimes of lO^-lO'' years (Equation [7|), kinetic energy is very efficiently converted 
to radiation, compared to nuclei with interaction timescales of 0(10^ y) (see 
Section [23]) . In addition, after 0(10 kyr), decelerating supernova shocks can no 
longer confine the highest energy particles, cutting off emission at VHE energies, 
whereas a powerful pulsar may drive a PWN significantly longer. 

Since electrons suffer energy-dependent energy losses as they convect /diffuse 
away from the pulsar, both synchrotron X-ray and IC 7-ray spectra should evolve 
with increasing distance from the pulsar. Indeed, this is observed in X-rays for a 
number of PWN, and in a single object at TeV energies. 

Table [2] lists 7-ray sources which are almost certainly associated with PWNe. 
Criteria for identification as a PWN include positional coincidence with a known 
pulsar powerful enough to potentially power the source, and the detection of an 
X-ray PWN. In some cases a radio SNR shell is also present providing a potential 
contribution to the 7-ray emission. Figure [9] shows some of the best 7-ray PWN 
candidates, with characteristics which were surprising at the time of discovery: 
(a) the large size of the 7-ray sources — typically ~ 0.5°, corresponding to a few 
tens of parsecs for typical few kpc distances, sometimes an order of magnitude 
larger then the corresponding X-ray PWNe and (b) the displacement of 7-ray 
sources from the pulsars by as much as the source radius, frequently putting the 
pulsar at the edge of the nebula. That the association is nevertheless significant 
can be demonstrated in two ways: a statistical study of sources in the H.E.S.S. 
Galactic Plane Surve y shows a significant 7-ray excess near powerful pulsars 
( Carrigan et al. 20071 ) an d in one case — HESS J 1825— 137 — energy-dependent 
morphology is observed ( Aharonian et al. 2006el ). with the source shrinking to- 
wards the pulsar for higher-energy 7-rays, as expected due to radiative cooling 
of electrons convecting away from the pulsar. 

The difference in size between 7-ray PWNe and X-ray PWNe can be attributed 
to the difference in energy of the electrons responsible for the radiation. For in- 
ferred B-fields at the PWN core of some 10 /xG, electrons of many tens of TeV 
are required to produce X-rays in the keV range (Equation [TT]). For these fields, 
cooling timescales are of order kyr, resulting in a modest range of particles. Elec- 
trons emerging from this inner region of the PWN still have energies sufficient to 
produce TeV 7-rays in interactions with IR and optical photons in the KN regime 
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(Equation [T0|) , but in the assumed /uG fields in the outer regions of the PWN 
the synchrotron peak is shifted into the optical, where it cannot be detected for 
all practical purposes. This also implies that VHE 7-ray emitting electrons, with 
characteristic cooling times of some 10 kyr, are accumulated over a much longer 
history of the pulsar, for medium-aged pulsars reaching back to the birth of the 
pulsar, where the energy output rate was one to two orders of magnitude higher. 
As a result, the 7-ray luminosity can reach and even exceed the current spin-down 
luminosity of the pulsar; an equilibrium is only reached once the pulsar age is 
large compared to cooling timescales, as is typically the case for X-ray emitting 
electrons. An exai nple for a resulti ng SED is shown in Figure [TU] for the object 
HESSJ1640-465 (|Funk et al. 2007l l. In the model shown here, electrons, up to 
20 kyr old are responsible for the detected VHE 7-rays, with their synchrotron ra- 
diation peaking (undetectably) in the visible. Young electrons, recently injected 
into the PWN, generate the compact X-ray nebula seen with XMM-Newton. 

Displacements between pulsars and their 7-ray PWNe appear to be fairly com- 
mon, but their origin is not fully understood; at least in the two cases where the 
proper motion of the pulsar is known (VelaX and HESS J1825— 137), the dis- 
placement between nebula and pulsar is almost orthogonal to the pulsar motion, 
eliminating the explanation that the pulsar was created with a significant kick, 
leaving a 'relic' PWN behind. The origin of the displacement is most likely a 
consequence of the environment, with dens ity gradients deforming the evolutio n 
of the SNR shell and hence also the PWN (jBlondin. Chevalier fc: Frierson 200ll ). 
Alternatively, target photon fields for IC scattering may be enhanced locally due 
to stars or star clusters. 

The sample of 7-ray PWNe is not large enough to systematically assess how 
PWN properties depend on pulsar properties. However, among the PWNe de- 
tected so far, two trends seem to emerge, with some caveats concerning selection 
bias (see e.g. Mattana et al. 2009l l: (a) the ratio of 7-ray luminosity to spin-down 
loss tends to decrease with increasing spin-down loss E (or decreasing pulsar age) , 
and the fraction of energy radiated in X-rays increases, and (b) very energetic 
pulsars (such as the Crab) tend to have very compact nebulae. A plausible expla- 
nation is that for rough equipartition between particle energy and magnetic fields 
in the nebula, the fields increase with E, which implies that the ratio of X-ray 
to IC luminosities increases with E and that the lifetime and range of electrons 
decrease. Furthermore, due to the strong correlation between pulsar age and 
spin-down luminosity, high spin-down pulsars tend to be young, such that the 
population of slower-cooled 7-ray emitting electrons is still increasing whereas 
X-ray emitting electrons have already reached their equilibrium output. The size 
of evolved 7-ray nebulae tends to saturate at some 10 pc (Figure E]). This may be 
because the nebulae themselves are confined by the ambient medium, or because 
radiative losses and convection timescales are such that multi-TeV electrons die 
out. 

By far the best-studied PWN is the Crab Nebula (jHester 20081 ). Its broad- 
band SED exhibits overlapping synchrotron and IC spectra, the first extending 
from radio into the EGRET energy range, the second from EGRET energies to 
~100 TeV (Figure [TT]) . Recently, the MAGIC collaboration succeeded in de- 
tecting, for the first time, a pulsed component in the VHE energy range, above 
25 GeV, hence excluding a pure polar-cap scenario as the origin of the pulsed 
emission ( Aliu et al. 20081 ). 
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4.3 Compact object binary systems as 7-ray sources 

The physical environment inside a close binary system (or an eccentric binary 
close to periastron) is radically different to that of the diffuse ISM relevant to 
SNRs and PWN. This environment is characterized by very high radiation densi- 
ties 0(lergcm~^) of rather high frequency photons 0(1 eV) and high magnetic 
fields (mG-G). The consequence of this environment for relativistic electrons is 
that rapid cooling is inevitable. In the case that the radiation pressure dominates, 
the cooling of TeV electrons will occur in the KN regime, with implications for 
the spectral shape discussed in Section 12.21 

All relevant timescales in such a system are short in comparison to the length of 
observation programs (typically years), for example the acceleration and cooling 
time for relativistic electrons and the orbital period. It is therefore expected that 
7-ray emitting binaries will be variable point-like objects if electrons dominate 
the 7-ray production. Hadron accelerating binaries may produce steady and 
extended 7-ray emission if the protons and nuclei can escape the production 
region without significant energy losses. If 7-rays are produced inside the binary 
system, then the assumption of free escape is normally invalid. In the presence of 
intense radiation fields, 7-7 interactions produce e''~/e~ pairs which are likely to 
IC scatter, leading to electromagnetic cascades. This effect, combined with KN 
IC cooling can lead to 7- ray spectral shapes radically different from tho se seen 



in diffuse sources (see e.g. iKhangulvan. Aharonian &: Bosch-Ramon 2008 ). 



The most obvious energy source inside a binary system in which one member is 
a compact object (neutron star or black hole) is accretion. Particle acceleration 
in jets produced by accretion onto a compact object is well established in extra- 
galactic objects and in the early 90s Galactic analogues to AGN were discovered 
and dubbed Micro-quasars ( Mirabel &: Rodriguez 19941 ). The internal and exter- 



nal shocks associated with such jets provide potential sites for particle accelera- 
tion. The primary alternative power-source is the collision of stellar/neutron-star 
winds, with acceleration at the termination shock of the wind, in a compressed 
version of a PWN. Of the three well-established systems, one, PSRB1259— 63/SS 
2883, is a unambiguously a binary PWN. For the other two systems (LS 5039 
and LSI +61 303) both PWN and micro-quasar scenarios have been extensively 
discussed. These objects are briefly descri bed below. Note that a much more 
complete historical account can be found in lAharonian et al. r2008hl ). 



The most detailed TeV measurements so far exist in the case of LS 5039, a 
3.9 day period syste m of a ^20Mc^ (06. 5V type) star and a compact companion 
of mass 3.7t\iMQ (iGasares et al 200^ ). Figure [12] shows the flux and spectral 



index of the V HE emission of LS 5039 a s a function of orbital phase, as measured 
usmg H.E.S.S. (I Aharonian et al. 2006al V As the distance between the stars varies 



by a factor ~2 around the orbit, the observed modulation of flux and spectrum 
with period is not unexpected. The maximum flux occurs at inferior conjunc- 
tion, the point where 7-7 absorption is expected to be at a minimum. However, 
the observed modulation is certainly not consistent with the simple-minded ex- 
pectation for such absorption (and/or cascading). As discussed in Section [231 
the peak absorption should occur at 0.9/£^T,eV = 300 GeV (as the thermal pho- 
tons from the companion star have kT ~ 3 eV), whereas there is apparently no 
modulation of flux at this energy. It seems likely that effects such as the angu- 
lar dependence of the IC cross-section, adiabatic losses and the dependence of 
acceleration efficiency (and also the efficiency of injection into the acceleration 
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process) and maximum energy, on distance between the massive star and com- 
pact object (i.e. changing accretion rate, shock velocity...) must be taken into 
account, together with the geometry of the system, to reach an understanding 
of the physical processe s at work. As LS 5 039 appears to host a bi-polar radio 
jet with a speed 0.2c ( Paredes et al. 2o"ool ) the accretion powered micro-quasar 
scenario has been most extensively discussed, but despi te extensive theoretical 
work (see Khangulvan. Aharonian &: Bosch-Ramon 20081 . and references therein) 
the acceleration site and nature of the compact object are still unclear in this 
system. 

LSI -1-61303 is a longer period (26.5 day) system with a lower mass compan- 
ion (~10 Mq, type BOVe) which h as been extensively observed by the MAGIC 
and VERITAS collaborations (see lAlbert et al. 2008d . and references therein). 
As in the case of LS 5039 VLBI radio observations had revealed extended jet- 
like radio emission in LS I -|-61 303, lead ing to its classification as a micro-quasar . 
However, more recent VLBI data from Dhawan. Mioduszewski &: Rupen (20061 ) 
suggest that the radio structure rotates with orbital phase, as might be expected 
for the cometary emission of a shocked pulsar-wind encountering a (higher mo- 
mentum) stellar wind from the massive companion. 

In the remaining well-established 7-ray binary system the pulsar nature of 
the compact object is well established by radio pulsation measurements. PSR 
B1259— 63 and its B2e com panion SS 2883 form a highly eccentric binary with an 
orbital period of ~3.4 years (I Johnston et al. 19921'). TeV emission from the systein 
close to its periastron passag e was predicted by Kirk. Ball &: Skiaeraasen (19991 ) 
and observed using H.E.S.S. (|Aharonian et al. 2005bl ~ PSRB1259— 63 is power- 
ful and close enough {E ~ 10^^ erg s~^, d ~2.3 kpc) that TeV emission might 
have been expected for a classical PWN. The TeV detection around periastron is 
usually attributed to the boost in IC emission from the strong radiation field of 
SS2883. 

The detection using the MAGIC telescope of a single flare from Cyg X-1 is ex- 
tremely important as in this system there is no doubt about the nature of the com- 
pact object - it is a 21 lb 8Mq black hole. The ~80 minute fl are has a significanc e 
of approximately 4.1 a after accounting for statistical trials ( Albert et al. 2007dl ). 
At present there is therefore strong evidence for, rather than proof of, VHE emis- 
sion from Cyg X-1. 

Given the relatively deep survey of most of the Galactic volume using H.E.S.S., 
it is natural to ask if additional candidates exist for 7-ray binary systems. Of 
the ~20 unidentified VHE sources, very few are point-like in nature, as would 
be expected both by analogy with identified s ources. By far the best candidate 
for a new 7-ray binary is HESS J0632-b057 (I Aharonian et al. 2007 j ). a point- 
like TeV source coincident with a Be star and an EGRET source. Follow-up 
observations of this object with XMM- Newton led to the d iscovery of a variable 
X-ray source coincident with the star ( Hinton et al. 20^ ). If all these objects 
are truly associated then the SED resembles that of the known 7-ray binaries. 
This discovery suggests that a population of 7-ray binaries exists with somewhat 
lower X-ray, radio and 7-ray fluxes then the 3 well established systems, which 
have hence so far escaped detection. 
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4.4 Stellar clusters and stellar winds 

All known Galactic sources of VHE 7-rays are associated (directly or indirectly) 
to massive star formation. Both the end-points of the massive stellar lifecycles, 
SNRs and pulsars, and high mass stars with compact companions (HMXBs) are 
TeV emitters. It is natural to consider whether massive stars can accelerate par- 
ticles to TeV energies in the absence of a compact object. The idea of particle 
acceleration at the shock front formed by colliding stellar winds in a binary system 
of two massive stars has been developed over the last five years (see for example 
Benaglia &: Romero 20031 : Pittard &: Doughertv 2006 ; Reimer. Pohl &: Reimer 20061 ) .| 



The dis covery of non-thermal hard X- ray emission from the massive binary Eta 
Carina ( Levder. Walter &: Rauw 20081 ) has strengthened the case for high-energy 



particle acceleration in these systems. 

So far there are no unambiguous VHE detections of individual colliding wind 
systems. However, extended TeV 7-ray emission has been detected from in and 
around Wes ter lund 2, the second la rgest concentration of massive young stars in 
our galaxy ( Aharonian et al. 2007d ). It seems plausible that this source is pow- 



ered by the collective effect of stellar winds within the cluster. The association of 
massive stars Cyg 0B2 has also been suggested as the co unterpart of the uniden - 
tified Galactic plane source discovered using HEGRA (|Aharonian et al. 20021 ). 



However, as essentially all known types of Galactic 7-ray source are associated 
(directly or indirectly) with high-mass star formation it is plausible that these 
emission regions are associated with a single object within the cluster, such as a 
PWN, which has not yet been identified at other wavelengths, rather than the 
cluster whole. 



4.5 Unidentified VHE 7-ray sources 

Roughly one third of the ~60 Galactic TeV sources have no compelling counter- 
part at other wavelengths. In several cases sensitive follow-up X-ray and radio 
observations have failed to identify these sources. Due to the apparent lack of 
emission of these TeV sources at lower frequencies, they are sometimes referred 
to as "dark accelerators" . The main questions which arise for this population are 
1) do they represent a new class of objects, or are they members of the known 
classes? and 2) is the emission hadronic or leptonic in origin? The identified 
VHE sources have low-frequency counterparts with non-thermal emission usually 
attributed to the synchrotron process. The lack of synchrotron counterparts to 
the unidentified sources may be taken as a suggestion that the 7-ray emission is 
produced by hadrons rather than leptons in most of these unidentified sources. 
However, there are several complications to this simple picture. 

The distribution of sizes and spectra of the unidentified sources are rather sim- 
ilar to those of the identified sources. However, these properties may be common 
to most TeV sources on rather general physical grounds. For example, the rela- 
tively fast diffusion and slow energy losses of >TeV particles make TeV sources, 
in general, rather large. The expected size of an old (i.e. in equilibrium) source of 
electrons cooled by synchrotron emission and with Bohm diffusion is rpc ~ SOS^q 
(from equations [1] and [7|) . Typical ISM magnetic fields of a few fiG therefore lead 
to sources with a scale of a few parsecs, with apparent radii of ~0.1° if located at 
typical distances of a few kpc. Unfortunately, the difficulty of identifying a source 
increases rapidly with its angular extent (unless closely matched in morphology to 
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a bright object in another waveband as in the case of the 7-ray SNR shehs). The 
sensitivity of existing X-ray and radio surveys for very extended ~0.5° objects is 
Hmited, and along the Galactic plane several candidates typically exist per degree 
of longitude. Source confusion can therefore be a ma jor difficulty. Unidentified 



sources such as HESS J1745— 303 (jAharonian et al. 2008c ) lie in regions with a 



high density of candidates and very likely have unresolvable contributions from 
several objects. Even with better angular resolution the inherently diffuse nature 
of CR sources may make such regions very difficult to disentangle. 

Leptonic scenarios for unidentified TeV sources are often dismissed as requiring 
very low magnetic fields to avoid bright synchrotron counterparts. However, 
even for 0(10 fiG) magnetic fields it is plausible that the synchrotron emission 
accompanying the VHE IC emission is largely confined to the FIR~UV band, 
and hence almost impossible to detect for even the brightest sources (e.g. i/Fv 
~ ]^Q-io gj.g cm~^ s~^) if they are extended on typical ~0.2° scales. X-ray 
synchrotron emission with such B-fields requires higher electron energies than 
needed to produce the ~ 1 TeV IC 7-rays. Furthermore, it is also possible to 
avoid radio counterparts if there is a low-energy cut-off in the injection spectrum 
of electrons, quite plausible in the case of, for example, pulsar wind nebulae. 

The lack of a radio counterpart may be harder to explain in the context of 
a hadronic scenario. Secondary electrons down to GeV energies are inevitably 
produced in the same p-p collisions that produce ir^ 7-rays. The accumulation of 
these electrons (assuming an E""^ injection spectrum) will lead to a synchrotron 
energy flux 

uFsync ~ 3 X 10-\uGH.)^/\B/WfiGf/^{t/10'yr)uFT:,y (17) 

at intermediate energies where the low energy turn-over in the electron distri- 
bution can be ignored and cooling is not important. In the saturated regime at 
high frequencies (tcool age) the synchrotron energy flux is directly proportional 
to the 7-ray flux and Fsync ~ O.lSF^o (see Figure [3]). For a typical TeV source 
with an energy flux of ~ 10~^^ ergcm"^ s~^, the resulting extended X-ray source 
would be difficult to detect with current instruments. 

In this context it is useful to consider some examples, both of sources which re- 
main unidentified and those whi ch were initially unident ified and now have likely 
counterparts. HESS J1303— 631 ( Aharonian et al. 2005d ) is an example of a TeV 



source which originally appeared to be without any compelling candidate at lower 
frequencies. At the time of discovery a PWN associated with PSR J1301— 6305 
was considered rather unlikely due to the required extremely efficient conversion 
of spin-down power to TeV emission and the large offset and rms size of the TeV 
source (both ?a 0.16°). In the light of the TeV PWN discoveries discussed in 
Section [121 PSR J1301— 6305 is now considered a compelling candidate. Given 
an age of ~ 10^ years, evolutionary effects are likely to enhance the TeV signal 
(see Figure [TUj) . Furthermore, the pulsar is quite plausibly much closer than its 
nominal dispersion measure distance of 15 kpc, making the offset and size very 
typical for a TeV PWN. 

The second class of "no-longer unidentified" objects contains those where follow-| 
up observations have led to the discovery of new pulsars, PWN and/or SNRs. 
HESS 1813— 178 is an example of such an object, rapidly established as a new 
composite SNR following its discovery, based on new and archiva l data from the 



VLA, INTEGRAL, XMM and Chandra (see iHelfand et al. 20071 . and references 



therein). HESS 1813— 178 is, however, an unusual case. It is a rather bright 
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(w 3 X lO^"*^^ erg cm~^ s^'"'^), compact (w 2') TeV source with a bright X-ray 
counterpart (AX J1813— 178). Such rapid assignment to an existing source class 
seems unhkely for the remaining unidentified sources. Figure [13] shows four exam- 
ple VHE sources for which there are no plausible candidates at lower frequencies. 
All lie within half a degree of the Galactic plane and all are significantly extended 
beyond the instrumental PSF. 



4.6 Active galaxies 

The galaxy Mrk 421 was the second VHE 7-ray source detected (IPunch et al. 1992l l 



the number of extragalacti c VHE sources has n ow risen to well over 20, with red- 
shifts up to 0.536 (3C 279, 1 Albert et al. (;2008bl )^. In order of increasing redshift. 



the VHE 7-ray emitters include M87 (z = 0.004), Mrk 421 {z = 0.030), Mrk 501 
{z = 0.034), lES 2344+514 {z = 0.044), Mrk 180 {z = 0.045), lES 1959+650 
{z = 0.047), PKS 0548-322 {z = 0.069), BL Lacertae (z = 0.069), PK S 2005-489 
{z = 0.071), RGB J0152+017 ( z = 0.08) (jAharonian et al. 2008bl ). W Comae 



{z = 0.102) (jAcciari et al. 20081 ). P KS 2155-304 (z = 0.1 16), H 1426+428 



0.129), lES 0806+524 {z = 0.138) (IV. Acciari et al. 2009 ). lES 0229+200 {z 



0.139), H 2356-309 {z = 0.165), lES 1218+304 {z = 0.182), lES 0347-121 
{z = 0.188), lES 1101-232 (z = 0.186), lES 1011+469 (z = 0.212), 3C 279 
{z = 0.536), PG 1553+113 {z unknown). A compilation of emis sion characteris - 
tics and remaining references for these objects can be found in IWagner f2008l ). 



All these objects harbor Active Galactic Nuclei (AGN), where a supermassive 
black hole with a mass from millions to billions of solar masses accretes mat- 
ter and powers jets — collimated highly relativistic outflows. Unlike Galactic 
VHE sources, all objects appear point-like given the ~5' resolution of lACTs, 
and all VHE source positions are consistent with the nominal location of th e 
AGN. With the exception of the radio galaxy M87 ( Aharonian et al. 2006fl ). 



all belong to the blazar class where a jet points towards the observer. High- 
frequency peaked BLLac objects (HBLs) dominate the sample. The few ex- 
ceptions are the flat spectrum radio quasar (FSRQ) 3C 279, the intermediate- 
frequency peaked BLLac object (IBL) W Comae and the low-frequency peaked 
BL L ac object (LBL) B LLacertae itself. Some obj ects — most notably Mr k 421 
(e.g. iFossati et al. 20081 ') and PKS 2155-304 (e.g. lAharonian et al. 2007^ ^1. see 



Figure O — exhibit burst-like variability on timescales of a few minutes to a 
few tens of minutes; emission from AGN such as M 87 or 3C 279 appears vari- 
able on day timescales. Despite the vastly larger distances, some AGN, at the 
peak of flares, outshine the strongest galactic sources by more than a factor 
of 10. While low-statistics measurements of VHE blazar spectra are consis- 
tent with power-laws, well-measured spectra are generall y significantly curved , 
steepening with incr easing energy, for exani ple Mrk 421 ( Krennrich et al. 200 ll ) 



and PKS 2155— 304 ( Aharonian et al. 2007bl ). Spectral indices tend to increase 



with source distance, at least partly due to absorption of high-energy 7-rays 
on infrared intergalactic photon fields, see Section 14. 7t but perhaps also related 
to the fact that distant AGN must be intrinsically brighter to be detectable. 
Spectral indices vary with flux for some sources, with a tendency for sp ectra 
to harden with increasing activity (see for example iKrennrich et al. 20021 ) . All 



VHE-detected blazars are also relatively strong radio and X-ray sources; in fact, 
observ ation targets are typically selec ted on the basis of their output in these 
bands ( Costamante Ghisellini 20021 ). Optical measurements have also proven 
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useful in tagging hig h activity states of A GN for targeted observations by VHE 
instruments (see e.g. Albert et al. 2006al ). The broad-band spectral energy dis- 
tributions of AGN exhibit a double-humped shape, with a high-frequency peak 
in the GeV to TeV regime, and a low-frequency peak in the optical to X-ray 
regime (Figure \T5\i . Simultaneous me asurements in the X-ray and VHE 7-rays 
bands reveal — with a few exceptions ( Daniel et al. 20051) — strong correlations 
between the X-ray and 7-ray fluxes (e.g. Fossati et al. 20081 ). suggesting a com- 
mon electron population as the origin of the radiation, with synchrotron X-ray 
emission and inverse Compton 7-ray emissior i. Similarly, av erage X-ray and 7-ray 
luminosities show a pronounced correlation ( Wagner 20081 ). 

In these blazars, VHE 7-ray emission is thought to arise inside the jets; mod- 
els typically assume a spherical 'blob' of high-energy particles moving along the 
jet axis, with flares created when high-energy electro ns are freshly injected into 
the blob (e.g. Tavecchio. Maraschi &: Ghisellini 19981 ). The relativistic motion of 
the blob beams and boosts the emission and reduces power requirements by a 
factor (5^, with 5 denoting the Doppler factor 6 = (r(l — l3cos9))~^ where F is 
the Lorentz factor describing bulk motion of the jet and 6 is the angle between 
jet and observer. The blob size is usually estimated by causality arguments from 
the minimum timescale tyar of variability, R ~ ct^ar^ ■ Minimum Lorentz factors 
are then obtained from the requirement that the energy density in the blob rest- 
fram e must be low enough to make the blob reasonably transparent for 7-rays 
(e.g. iDondi fc Ghisellini 19951 ). and are in the range of 10 to almost 100. Mag- 
netic flelds in the blob, as well as parameters of the electron spectrum, can then 
be determined from the relative spacing and height of the two peaks of the SED 
and from the requirement that the cooling time of electrons matches the flare 
timescale, and are typically in the range 0.01 G to 1 G. For the shape of the elec- 
tron spectrum, a broken power-law is frequently employed. Parameter estimates 
should, however, be taken with a grain of salt; well-sampled AGN light curves are 
com posed of a range of Fouri er components with a red-noise spectral distribution 
(e.g. Aharonian et al. 2007bl ). without a highest frequency which could be unam- 
biguously as sociated with the size of a blob. In Sync hrotron-Self-Compton (SSC) 
models (e.g. Tavecchio. Maraschi fc: Ghisellini 19981 ) synchrotron photons gener- 
ated in the blob provide the dominant target for IC up-scattering by electrons, 
resulting to flrst approximation in a quadratic increase of 7-ray flux with electron 
number and X-ray flux. Detailed modelling of broad-band spectra, however, also 
indicates the relevance of external target radiation components from other regions 
of the jet and from the high-te mperature gas surrounding the central engine (e.g. 



Katarzvhski. Sol fc Kus 200lh (Figure [E]) . 



Phenomenological models frequently leave open how particles in the blob are 
accelerated; possibilities include shock-wave acceleration in MHD turbulence in 
the jet, centrifugal acceleration of particles along rotating rn agnetic fleld lines 
near the base of the jet, or shear ( Rieger &: Mannheim 20o3 ). However, X-ray 
observations show that acceleration of high-energy particles must occur along the 
jet, since travel times from the base of the jet to X-ray emitting knots exceed 
cooling times. Many fundamental aspects of AGN jets and of particle accelera- 
tion in these jets are poorly understood, including the mechanisms which launch 
the jets and their composition. Not even the character of the emitting particle 
pop ulations is flrmly established: cascades induced by high-energy protons (see 
e.g. Mannheim 19931 ) can successfully reproduce most features of AGN SEDs, 
with the possible exception of very fast variability, which is hard to model due 
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to the long energy loss timescales of protons in comparison to electrons. VHE 
measurements are a promising method for probing jet properties and ultimately 
studying the environment of supermassive black holes. 

The discovery of VHE 7-ray emission from M 87 established radio galaxies as a 
new class of source where 7-rays are emitted at significant angles to the jet, not 
relying on Doppler boosting (and hence only observable f or nearby objects). The 



observation of fast - day scale - variability of the emission (jAharonian et al. 2006: 



Albert et al. 2008al ) excludes the extended jet or the radio lobes as sources and 



imply a compact emission region with a size comparable to the radius of the 
SMBH at the center of M 87, most likely the nucleus itself or possibly the knot 
HST-1 in the inner jet. 

4.7 Probing background radiation fields with 7-rays 

The energy spectra of extragalactic 7-ray sources are n iodified by interactions 
with the diffuse extragalactic background light (EBL, see lHauser &: Dwek f 200 il l 



for a review) , see Section 12.51 and Figure HI Since the absorption cross section 
peaks near threshold, there is an approximate mapping between the wavelength of 
an absorbing EBL photon and the energy of an interacting 7-ray, -BxeV ~ 0.7A^m- 
The level of the EBL, in particular at the relevant IR wavelengths, is very dif- 
ficult to determine by direct observation, due to overwhelming foreground light 
sources, hindering the interpretation of AGN spectra due to the uncertain cor- 
rection for EBL absorption. On the other hand, the EBL is important in its own 
right and the information provided by this TeV absorption is potentially very 
useful. The EBL represents the integrated, red-shifted, emission from all epochs 
of the evolution of the Universe. The EBL, with its peaks around l^m from 
starlight and around 100//m due to starlight reprocessed by dust (Figure [T6|) . has 
embedded within it information on the history of galaxy and star formation in the 
Universe. Assuming 'plausible' shapes for blazar spectra, one can use absorption 
featur es in 7-ray spectra to constrain th e level and spectral distribution of the 



EBL (jStecker. de Jaeer fc Salamon 19921 ). The 'plausible' assumptions usually 



include a 7-ray spectral index F^ > 1.5 at the source, and emitted spectra which 
are power-laws, possibly curved with the index increasing — but not decreasing 
— with increasing energy. These assumptions are supported by nearby sources, 
which are less affected by EBL absorption. However, it can not be firmly excluded 
that the more distant sources used to 'meas ure' the EBL are subject to a selectio n 
bias and differ in their intrinsic spectra (e.g. Stecker. Baring fc Summerlin 200?! ). 



EBL de-absorption with an optical depth t{E) transforms a measured 7-ray fiux 
^oiE) into an intrinsic fiux at the source ^i{E) = ^o{E)e'^^^\ A rapid increase 
of t{E) at higher energies, as predicted in some EBL scenarios, will decrease 
the intrinsic spectral index and may confiict with the assumptions concerning in- 
trinsic spectra. Sources with higher redshifts and correspondingly large t{E) and 
with hard measured spectra and wide energy coverage will provide the most strin- 
gent constraints. On the basis of spectra measured for lES 1101—232 {z = 0.186) 
and H 2356— 309 {z = 0.165), upper limits for the EBL density at optical/near- 
IR wavelengths were derived which were within a factor 1.5 of the lower limits 
provided by the int egrated light of resolved galaxies, and below the EBL level 
previously assumed ( Aharonian et al. 2006cl ). Strong constraints on the shape of 



the EBL spectrum in the 5 -10 //m region were obta, ined based on the spectrum of 
lES 0229+200 {z = 0.139) (jAharonian et al. 2007al ). The dection of VHE 7-rays 
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from 3C 279 {z = 0.536) ( Albert et al. 2008bl ) provided another demonstration 



of low EBL levels and hence the relatively high transparency of the Universe to 
7-rays. Limits derived for the EBL density often include — beyond the assump- 
tions on intrinsic source spec tra — certain minima, l assumptions about the shape 
of the EBL spectrum itself. Mazin &: Raue (200?! ) have developed an approach 



where arbitrary EBL shapes are tested and constrained, confirming earlier EBL 
limits (Figure [T6]l . EBL determinations using absorption of VHE 7-rays have 
driven development of models for EBL fo rmation, and the latest models (see e.g. 
Franceschini. Rodighiero fc Vaccari 20081 ) and calculations, based on hierarchi- 



cal structure formation and employing detailed models for galaxy formation and 
evolution and the reprocessing of starlight by dust, repro duce the low EBL levels 
measured (see e.g. Primack. Gilmore &: Somerville 20081 ). 



4.8 Other extragalactic source classes 

A range of extragalactic objects beyond AGN is expected to emit VHE 7-rays at 
some level, such as normal galaxies, starburst galaxies, galaxy clusters, GRBs, 
and the sources of ultra- high-energy cosmic rays (UHECR); most of these are 
predicted to emit at a level not far below the sensitivity of current instruments, 
and many of these classes raise interesting new science topics. 

Nearby galaxies. Diffuse 7-rays from nearby normal galaxies are produced 
in CR interactions with the interstellar medium; their detection - combined with 
estimates of supernova rates in these galaxies - would allow concepts of CR accel- 
er ation and propagation in galaxies to be tested. Expected fluxes are discussed 



m 



Pavlidou &: Fields (20011 ). applied to the EGRET energy range. However, 



detection of local group Galaxies at VHE energies will be difficult even for next- 
generation instruments. 

Starburst galaxies and ultraluminous infrared galaxies. Starburst galax-| 
ies with strongly enhanced supernova rates and enhanced gas density, and their 
more extreme cousins the ultralumious infrared galaxies, should provide VHE 7- 

ray fluxes not far below current sensitivities, providing anot her test of concepts of 

CR accele ration. Amongst the best candi dates are NGC 253 ( Domingo-Santamaria fc Torres 2005l )| 



and M 82 (iPersic. Rephaeh fc Arieh 2008i l: i n fact, NGC 253 w as claimed as a TeV 



emitter by the CANGAROO collaboration (jltoh et al. 20031 ). but this detection 



was not c onfirmed by later , more sensitive, observations with H.E.S.S. and CAN- 
GAROO (jltoh et al. 20071 ). Bounds on emission from the ultralumious infrared 



galaxy Arp 220 have been reported by the MAGIC collaboration ( Albert et al. 2007bl ).| 



Galaxy clusters. Clusters of galaxies are the largest gravitationally bound 
systems in the Universe, and should contain a significant non-thermal particle 
population, fed by particle acceleration in accretion shocks during the assembly of 

the cluster, by supe rnova activity in cluster galaxie s , and by particle acceleration 

by cluster AGN (e.g. lBlasi. Gabici k Brunetti 2007l : IVolk. Aharonian k Breitschwerdt 1996l ).| 



The time required for CRs to diffuse out of a cluster is generally supposed to be 
larger than the age of the Universe. As a consequence the CR abundance in 
clusters provides a measure of the time-integrated CR production rate, and the 
CR spectrum is not softened by escape as it is, for example, inside our own 
galaxy. VHE 7-ray flux limits have been derived using t he Whipple telescope for 
the Perseus and Abell 2029 clusters ( Perkins et al. 200^ ). limiting the CR energy 



density to as little as 8% of the thermal energy density in the case of the Perseus 
cluster, depending on assumptions concerning the distribution of CRs relative to 
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gas in the cluster. Such a value is well within the range estimated for different 
acceleration mechanisms, indicating that detections may be possible in the near 
future. 

GRBs. Gamma ray bursts - GRBs - are usuall y explained by fireball models, 
with emission produced by relativistic shocks (e.g. Meszaros 20061 ). Both prompt 
and delayed 7-ray ernission at GeV and Te V energies has been predicted (e.g. 
Asano &: Inoue 20071 . Pe'er &: Waxman 20051 ). due to both leptonic and hadronic 



mechanisms. No VHE 7-rays have been detected from GRBs so far. Truly 
simultaneous observati ons are primarily p ossible with wide-field instruments such 
as MILAGRO (see e.g. lAbdo et al. 2007bl . for limits on 17 GRBs), at the expense 
of sensitivity and energy threshold in comparison to current Cherenkov telescopes. 
In one case — GRB 060602B — an (unusually soft) GRB occurred in the field of 
view of ongoing observatioi is of the H.E.S.S. teles cope, however, without giving 
rise to a detectable signal ( Aharonian et al. 20091 ). Upper limits were reported , 
for example, from the MAGIC collaboration for 9 GRBs (jAlbert et al. 2007cl ). 
with observations starting as early as 40 s after the burst, due to the fast slewing 
capabilities of the MAGIC telescope. A strong and nearby GRB, observed rather 
promptly, is needed to challenge GRB models. None of the limits obtained so far 
fall in this class and it looks as if patience will be required to obtain a VHE 7-ray 
detection of a GRB, even if a TeV emission component is present, due to the 
limited duty cycle and redshift reach of the sensitive instruments. An in-depth 
discussi on and detailed referen ces concerning high-energy 7-rays from GRBs is 
given in lBucklev et al. ^2008^ 1. 



Sources of UHECR. Nearby accelerators of ultra high energy cosmic rays 
(UHECRs) - within the sslOOMpc GZK radius - should also be sources of VHE 
7-rays. Firstly, interactions during the propagation of UHECR beyond the GZK 
cut-off give rise to cascades fee ding energy down to the TeV range (see e.g. 
Ferrigno. Blasi &: de Marco 20051 ). Secondly, VHE photon prod uction during ac - 



celeration of UHECRs is expected in some scenarios (see e.g. Levinson 2000l ). 
Again, fluxes for sources from which a few UHECR are detected are typically 
expected to lie at the lower edge of the sensitivity of current instruments. 



5 VHE 7-rays and astroparticle physics 

The areas of VHE 7-ray astronomy discussed so far concern primarily issues in 
astronomy and astrophysics. The development of VHE 7-ray instrumentation has 
also been driven to a significant extent by questions in the field of astroparticle 
physics, above all the indirect search for dark matter. Results obtained by VHE 
instruments also have relevance to fundamental physics, for example in tests 
for an energy dependence of the speed of light, as predicted in some models of 
quantum gravity. The former aspect will be addressed in more detail, the latter 
rather briefly. 



5.1 Searching for dark matter 



In the standard cosmological model (e.g. lFrieman. Turner &: Huterer 20081 ). which§ 

has emerged over the last decade, only about 4% of the energy density in the 
Universe is in the form of baryons, about 25% is cold dark matter (DM) and 
the remainder is dark energy. Evidence for the gravitational effects of DM is ob- 
served on all scales ranging from galaxies, where flat rotation curves require that 



Teraelectronvolt Astronomy 



27 



visible matter is embedded in extended DM halos, to clusters of galaxies. DM is 
essential in explaining structure formation in the Universe, driving the formation 
of potential wells in which baryonic matter accumulates, forming galaxies. Struc- 
ture formation favors cold dark matter — i.e., DM in the form of non-relativistic 
particles — since rapidly streaming DM would have smoothed inhomogeneities in 
the matter distrib ution on small scales. Sorn e form of weakly interacting massive 
particle (WIMPs ( Steigman fc Turner 19851 )) is required to explain the existing 
observations. 

Creation and annihilation of DM particles was in equilibrium in the early Uni- 
verse, but under the rapid expansion of the Universe the WIMP component was 
rapidly frozen out, the small interaction cross-section no longer sustaining suffi- 
cient annihilation rates. The current DM abundance emerges naturally assuming 
typical weak interaction cross sections and WIMPs in the mass range of some tens 
of GeV to TeV. Candidates include super symmetric particles, where the lightest 
supersymmetric particle is stable due to (assumed) R-parity conservation, and 
Kaluza-Klein particles arising in theories with TeV-scale extra dimensions. As 
a result of gravitationally-driven structure formation, DM particles form halos 
in which galaxies are embedded, with a pronounced peaks at galactic centers. 
Simulations of structure formation predict a 1/r (NFW) density profile near the 
center of galaxies, with a significant number of sub-halos scattered elsewhere 
within galaxies. However, details of the merger history of halos/galaxies as well 
as the influence of central black-hole formation may significantly infiuence the 
distribution of DM, resulting in large uncertainties in the predicted densitie s. 

Identification of the nature of DM constituents (e.g. Hooper &: Baltz 20081 ) has 
been a major theme of particle and astroparticle physics in the last one to two 
decades. Possibilities for the detection of cosmological DM include challenging 
direct-detection experiments, where DM particles — typically assumed to have 
a local density of a round 0.3 GeV cm~^ — scatter off nuclei, resulting in keV- 
scale recoil energies (jCaitskell 20041 ). Annihilation of DM particles in the galactic 
vicinity of the Sun will enhance CR electron and positron fluxes and in particular 
the positron/electron ratio, as well as the antiproton/proton ratio, depending on 
the specific annihilation channels. Finally — and most relevant for VHE 7-ray 
astronomy — annihilation of DM particles will generate an enhanced 7-ray fiux 
from regions of high DM density, the annihilation rate reflecting two-particl e in- 
teractions, proportional to the square of the density (see lBucklev et al. 2008al . and 
references therein). 7-rays may be produced directly in the primary annihilation 
process, by the decay of hadrons created via annihilation into fermion-antifermion 
pairs or vector bosons, or by IC scattering of electrons created in the annihilation 
or in the decay chains of annihilation particles. 7-rays from DM annihilation are 
identifled by their spectral and directional signature. The expected flux is gen- 
erally written as a product of a "particle physics" factor and an "astrophysics" 
factor, 

*.-(^«,(^))(it/^^0 m 

where < cif > is the velocity- weighted annihilation cross section of WIMP parti- 
cles of mass M, D^{E) the spectral yield per annihilation, and p the density of DM 
particles, with the annihilation flux being proportional to the line-of-sight (LOS) 
integrated squared density. For all decay modes the spectrum D^{E) cuts off at 
the mass of the DM particles, which annihilate effectively at rest. The spectral 
details depend very much on the decay modes. Two-body decays into 7-rays, such 
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as 77 or Z^^ or Hj, resulting in a monoenergetic line signature, are possible via 
loop processes, but are strongly suppressed in most scenarios. Dominant decays 
into quark-antiquark pairs or boson pairs, followed by hadronization and hadron 
decays, generate a broad featureless spectrum with an SED hump at about 10% 
of the WIMP mass (see Figure [T7j) . Decays into tau leptons result in a harder 
7-ray continuum. In one class of models, two-body decays are helicity suppressed, 
boosting diagrams with internal brer nsstrahlung such as W~^W~'y and a 7 -ray 



enhancement near the kinematic limit (jBringmann. Bergstrom fc Edsio 20081 ) , al 



most equivalent to line emission given the finite energy resolution of VHE 7-ray 
detectors (Figure [T5|) . 

The expected flux of annihilation 7-rays is proportional to the LOS integrated 
squared DM density. Searches for signatures of DM annihilation in the VHE 
regime therefore concentrate on objects with spikes in DM concentration. Fa- 
vorable objects include the Galactic center, dwarf spheriodal galaxies such as 
Draco, Ursa Minor and Ursa Major, or Willman 1 ( Strigari et al. 20081 ) as well 



as microhalos and DM enhancements p redicted around intermediate-m ass black 
holes (IMBHs) in the Galactic vicinity ( Bertone. Zentner &: Silk 20051 ). globular 



clusters, and supermassive black holes in external galaxies. A DM signal would 
have to be identified by a combination of spectral and directional signatures, 
with nearby sources exhibiting a narrow spike towards the center of the source 
combined with an extended tail. 

In the VHE domain, despite extensive searches, no DM candidate source has 
been identified so far. In fact, predicted fluxes tend to be well below the sensitivity 
of current instruments, requiring "boost factors" of several orders of magnitude to 
make signals detectable, arising from enhancements in the DM density. Detection 
of a DM signal from the most obvious candidate source, the Galactic Center, is 
hampered by the presence of a strong astrophysical source of 7-rays, responsible 
for a power-law spectrum of 7-rays extending beyond 10 TeV, which may hide 
a faint low-energy annihilation signal . Limits on annihila tion fluxes have been 



obtained e.g for the Galactic Center (jAharonian et al. 2006h ) , the dwarf galax - 



ies Draco (lAlbert et al. 2008dl.lWood et al. 20081 ^. Ursa Minor (IWood et al. 20081) 



and Sa.it tkrius Uhaxonia.et ;i. 2008^). for 1ocd .rou n galaxi es ^Wood et ab 20 ^ 1 



and for the globular cluster M15 ( Wood et al. 2008l l. DM annihilation around 



IMBHs in the Galactic halo predicts detectable fluxes eve n in the absence of 
large additional boost factors ( Bertone. Zentner &: Silk 20051) and limits res trict- 



ing the model parameter space have been obtained ( Aharonian et al. 2008gl ). but 
the process relies on a specific astrophysical scenario which may not be realized 
in nature. Interestingly, DM annihilation fiuxes from IMBHs depend only weakly 
on the annihilation cross-section, since close to an IMBH the DM density is so 
high that annihilation limits the peak density, so that the resulting p^a varies 
only as cj^/^. 

5.2 Violation of Lorentz invariance and quantum gravity 

In the framework of quantum gravity models, Lorentz invariance may be vio- 
lated such that the propagation spee d of radiation depends on the photon energy 
and possibly on its helicity (see e.g. Amelino-Camelia et al. 19981 ). For energies 



small compared to the scale of quantum gravity, the photon speed c' can be 
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parametrized as 

c' = c( 1 + + ) (19) 




where Ep = 1.22 • 10^^ GeV is the Planck energy and and C sue free param- 
eters which are expected to be of order unity, unless specific symmetries forbid 
the corresponding terms; for example, there are classes of models with = 0. 
GRBs and flares from AGN can be used to perform "time of flight" studies, 
searching for an energy dependence of the peak time of arrival of the emission. 
The best sensitivity in terms of limits on ^ and C is obtained for short bursts 
with well-located peak emission, for distant objects and for large energy differ- 
ences. A positive detection of energy-dependent peak emission could of course 
also be attributed to mechanisms inside the source, where such effects can arise 
rather naturally, governed by acceleration and/or cooling time scales. To con- 
flrm that a detection implies violation of Lorentz invariance, one would need 
to demonstrate that the effect increases with source distance, and that the dis- 
tance dependence cannot be attributed to evolution of the source population. At 
VHE energies, the tec hnique was flrst a pplied by the Whipple group to a flare 
of Mrk421 {z w 0.03) (|Biller et al. 19991 '). resulting in a limit of |^| < 200. More 



recently, MAGIC reported an energy dependence in the peak position of a flare 
of Mrk501 {z ^ 0.03), finally quoting a limit |^| < 60 (jAlbert et al. 20081 ). In 



the giant fiare of PKS 2155— 306 {z ^ 0.12 ) detected using H.E.S. S., no energy 
dispersion was visible, resulting in |,^| < 17 ( Aharonian et al. 2008al ). VHE 7-ray 



observations are hence beginning to probe energy-scales up to a few % of the 
Planck energy. While such time-of-flight measurements provide the most model- 
independent tests of a possible energy dependence of the speed of light, it must be 
mentioned that alternative methods — which make some additional assumptions 
— are signiflcantly more sensitive, p roviding limits from |£ | < 10~^ for helicity- 
depe ndent speed modiflcatio ns (e.g. Fan. Wei Xu (200?! ) down to |^| < 10"^"^ 



(e.g. iGalaverni Sid (200^ ) based on decay kinematics. 



6 Conclusions 



VHE 7-ray astronomy can now legitimately claim to be a mature astronomical 
discipline, with resolved source morphologies and well-sampled light curves and 
energy spectra. Furthermore, as some objects have their peak energy output in 
the TeV range, VHE observations are clearly critical to our understanding of some 
astronomical objects. The primary success of the field so far is in the unambiguous 
identification of sites of cosmic particle acceleration. TeV photons are currently 
the only effective probe of relativistic hadrons in astrophysical environments. Two 
complementary techniques exist for ground-based 7-ray astronomy: lACTs and 
ground-based particle detectors. Both have significant further potential, with 
order of magnitude more sensitive instruments such as CTA and HAWC being 
planned. The lACT technique has the potential to cover the energy range from 10 
GeV to hundreds of TeV, with an angular resolution of better than 1' achievable 
at a few TeV. Together, the Fermi satellite and lACTs will soon provide unbroken 
sensitive coverage over 7 orders of magnitude in 7-ray energy. 

The number of VHE 7-ray sources is now over 80. More importantly, many 
different source types are represented. Supernova remnants and pulsar wind neb- 
ulae appear to dominate the Galactic population, but pulsed magnetospheric 
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emission from the Crab pulsar has now been detected from the ground and a 
number of binary systems have been estabhshed as TeV emitters. There is also 
TeV emission apparently associated with clusters of massive stars, perhaps as a 
consequence of stellar wind interactions. Variable emission from active galactic 
nuclei dominates the extragalactic TeV sky. A complex mixture of different mech- 
anisms appear to be at work in these TeV sources. Bulk flow appears to be the 
main energy source; rclativistic jets/winds in AGN/pulsars and non-rclativistic 
shocks driven by supernovae ejecta. Acceleration processes transfer energy from 
these bulk motions into populations of relativistic particles with power-law en- 
ergy distributions. The convection, diffusion and cooling of these particles then 
determine the properties of the high-energy radiation. For accelerated hadrons 
cooling does not usually modify the injection spectrum and the observed 7-ray 
spectral shape is close to that of the injected hadrons - possibly with modifica- 
tions due to energy dependent propagation. The situation for electrons is more 
complex, cooling can in general not be neglected and there is a balance between 
synchrotron and IC losses - complex spectral shapes can emerge, particularly in 
the case of IC dominated cooling. 

Galactic TeV sources are typically extended on 10 pc scales. Compact 7- 
ray sources require both compact accelerators and dense compact target mate- 
rial/radiation fields/magnetic fields, and such objects appear to be rare — with 
binary systems as the only established example. The cooling time for 10 TeV 
electrons in 10~^ G magnetic fields is ^ years - comparable to the active 
lifetime of the dominant sources, energetic pulsars and young SNR. Both bulk 
motion at 1000 kms~^ over this timescale, and Bohm diffusion of electrons in 
older (cooling limited) objects, naturally produce 10 pc scale TeV sources. For 
sources on this scale the detection sensitivity of lACTs is very competitive to 
that of X-ray telescopes for detection of synchrotron counterparts. 

The current extragalactic TeV sources are characterized by relativistic bulk 
motion and compact regions of acceleration and cooling, leading to variability on 
timescales down to minutes. For the blazar class, the boosting associated with 
this bulk motion allows us to see distant objects — with a current redshift record 
of 0.536. The combination of TeV data with observations of optical to X-ray 
synchrotron emission promises to be a powerful probe of the inner jets of these 
AGN. Extragalactic TeV sources also provide a useful tool to study the cosmic 
infra-red background and fundamental physics such as the potential violation of 
Lorentz invariance at high energies. 

The next decade of TeV astronomy should bring the source count to 0(1000) 
and increase the diversity of known Galactic and extragalactic sources. Predicted 
VHE fluxes for several classes of astrophysical object are close to detcctability 
with current instruments, for example clusters of galaxies, starbursts and GRBs. 
The indirect detection of Dark Matter is a major objective of the field and consid- 
erable phase-space is available for discovery via this channel. For known source 
classes more sensitive and precise observations will bring improved understand- 
ing, in particular a quantitative understanding of Galactic CR origin seems within 
reach of the next generation of I ACT detectors. To conclude, TeV astronomy has 
emerged as a powerful tool for high-energy astrophysics and looks set to become 
our primary window on the non-thermal universe in the years to come. 
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Instrument 


Lat. 


Long. 


Alt. 


Tels. 


Area 


Pixels 


FoV 


Thresh. 


Sens. 




n 


n 


(m) 




K) 




n 


(TeV) 


(% Crab) 


H.E.S.S. 


-23 


16 


1800 


4 


428 


960 


5 


0.1 


0.7 


VERITAS 


32 


-111 


1275 


4 


424 


499 


3.5 


0.1 


1 


MAGIC 


29 


18 


2225 


1 


234 


574 


3.5t 


0.06 


2 


CANGAROO 


-31 


137 


160 


3 


172 


427 


4 


0.4 


15 


Whipple 


32 


-111 


2300 


1 


75 


379 


2.3 


0.3 


15 


HEGRA 


29 


18 


2200 


5 


43 


271 


4.3 


0.5 


5 


CAT 


42 


2 


1650 


1 


17.8 


600 


4.8t 


0.25 


15 



Table 1: Properties of selected air-Cherenkov instruments, including two of histor- 
ical interest (HEGR A and CAT), t These instruments have pixels of two different 
sizes. Adapted from lffinton (20081 ). 



Object 


Type 


Method 


Flux 


Ref. 


PSR B1259-63 


Binary 


Pos/Var 


7* 


Aharonian et al. f2005b) 


LS 5039 


Binary 


Pos /Per 


3* 


Aharonian et al. f2006a) 


LSI-F61303 


Binary 


Pos / Var 


16* 


Albert et al. (2006b) 


RXJ1713.7-39046 


SNR 


Mor 


66 


Aharonian et al. (2004) 


Cassiopeia A 


SNR 


Pos 


3 


Aharonian et al. (2001) 


Vela Junior 


SNR 


Mor 


100 


Kataeiri et al. (2005) 


RCW86 


SNR 


Mor 


5-10? 


Aharonian & et al. (2008) 


SN 1006 


SNR 


Mor 


? 


Naumann-Godo & et al. (2006) 


Crab Nebula 


PWN 


Pos 


100 


Weekes et al. (1989) 


G0.9-F0.1 


PWN 


Pos 


2 


Aharonian et al. (2005d) 


MSH 15-52 


PWN 


Mor 


15 


Aharonian et al. (2005a) 


HESSJ1825-137 


PWN 


EDMor 


12 


Aharonian et al. (2006e) 


VelaX 


PWN 


Mor 


75 


Aharonian et al. (2006p;) 



Table 2: Selected galactic VHE 7-ray sources with well established multi- 
wavelength counterparts. Note that all these objects are X-ray sources. Fluxes 
are approximate percentages of the TeV flux from the Crab Nebula, * indicates 
variable emission. The final column lists the publications where a firm identifica- 
tion of the source was made. These associations were established through a range 
of methods, given here in abbreviated form: Pos: The centroid position of the 
VHE emission is established with sufficient precision that there is no ambiguity 
as to the counterpart. Mor: There is a match between the 7-ray morphology and 
that seen at other wavelengths. EDMor. Energy-dependent morphology is seen 
which approaches the morphology seen at other wavelengths at some limit, and is 
consistent with our physical understanding of the source. Var: 7-ray variability 
correlated with that in another waveband is observed. Per. periodicity in the 
7-ra y emission is se en, matching that seen at another wavelength. Table adapted 
from iHinton (20081 ) ■ 
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Figure 1: Top: The current catalog of VHE 7-ray sources plotted on the sky 
in Galactic coordinates. Bottom: the Milky Way vie wed in VHE 7-rays: the 
H.E.S.S. survey of the Galactic plane (reproduced from lHoppe fc et al. 20071 ). 
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10"' 10"^ 10"' 10 10^ 10^ 10' 10' 10" io'^ io'' 




Energy (eV) 

Figure 2: a) SEDs for radiation of mono-energetic 1/100 TeV electrons (red/blue 
curves): Synchrotron and IC (solid curves) and Bremsstrahlung (dashed curves). 
Three IC curves are shown for each primary energy: (from low to high) on the 
CMB {kT = 2.35 x 10"^ eV, 5 4 x 10-3/0-4), on dust-emitted FIR {kT = 0.02 
eV, b ^ 0.3/30), and on visible (star) light {kT = 1.5 eV, b ^ 20/2000). Note that 
for 100 TeV electrons scattering on optical photons the IC energy distribution is 
effectively a delta-function at 100 TeV. The curve normalizations are appropriate 
for a total particle energy of lO'*^ erg at 1 kpc distance in a magnetic field of 3 fiG, 
a matter density of 100 hydrogen atoms cm^^ and radiation fields of density 0.26 
eV cm~3 (CMB and FIR) and 1 eV cm"'^ (starlight), b) SEDs for 7-rays and 
synchrotron radiation of secondary electrons from strong interactions of mono- 
energetic protons. The magnetic field is increased to 30 to illustrate the effects 
of cooling and steady injection over 10^ yr (dashed curves 10^ yr) is assumed. 
The input energy is 10^^ erg. c) and d) - as for a) and b) but for cut-off power- 
law distributions of particles: dN/dE oc E~'^ exp —E /Ec with Ec = I TeV (red) 
and 100 TeV (blue). 
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C/D 




Energy (eV) 

Figure 3: SEDs for continuous injection and cooling of a population of electrons 
with an E~'^ injection spectrum and an exponential cut-off at 100 TeV. Solid, 
short-dashed and long-dashed curves show injection timescales of 10^ yr, 3 x 10^ 
yr and 10^ yr, respectively. The blue curves show synchrotron and IC emission 
in the case of synchrotron-dominated cooling, with B = 30/iG and the CMBR as 
target for IC. The red curves illustrate the effects of IC-dominated cooling with a 
lower magnetic field {B = 3fiG) and a higher energy radiation field (black-body 
photons with kT=1.5 eV, with density 1000 eV cm~^) where KN effects become 
important. 
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Figure 4: Optical depth (r) of 7-rays as a function of energy, for FIR (/cT = 
0.008 eV) and visible {kT = 0.8 eV) target photon fields, of column density 3 • 
10^^ eV cm~^ for FIR (1 eV cm^'^ over 10 kpc) and 3- 10^"^ eV cm^^ for visible (1 eV 
cm~^ over 1 Mpc). The effect of absorption on a non-thermal photon distribution 
(with dN/dE oc E'"^) is shown for comparison. Note that the absorption is 
constant for an E^"^ photon field. The green curves show the optical depth for 
pair production on the EBL for redshift s of 0.01, 0.1 and 1 for the model of 
Franceschini. Rodighiero &: Vaccari (20081 ). 
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Figure 5: A sketch of the imaging atmospheric Cherenkov technique showing 
the formation of an electromagnetic cascade for a 300 GeV photon primary, the 
production of Cherenkov hght, and the formation of an image in the camera of a 
Cherenkov telescope. Cherenkov light production for a proton initiated cascade 
is shown for comparison. Shower images produced by Konrad Bernlohr. 
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Figure 6: Four SNRs imaged in (dominantly) non-thermal X-rays (left) 
and resolved in VHE 7-rays with H.E.S.S. (right). a) RX J1713.7-394 6 
with 1-3 keV data from ASCA ( Uchivama. Takahashi &: Ah aronia n 200^ ). 
b) RXJ0852.0-4622 with EDSAT fl. 3-2.4 keV) (lAschenbach 1998ly cl 
RCW86 with 2-4 keV data from XMM-Newton (jVink et al. 2nnfil ^ d) 
SN1006 with C handra archive data ( 0-5-10 k eV). The H.E.S.S. dat a 
are taken from Aharonian et al. r2006bl . l2007dl l. lAharonian fc et al. 
Naumann-Godo fc et al. f2006l ). The white scale bars are 0.5° long. 
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Figure 7: The SED of RX J1713.7— 3946 at 7-ray ene rgies. Three curves are 
shown are shown in comparison to the H.E.S.S. data ( Aharonian et al. 2007el ) 
and EGRET upper hmit: a) the best-fit 7-ray spectrum arising from inter- 
acting protons wit h an energy distribution following a power-law with expo- 



nential cut-off, s ee iKelner. Aharonian Sz Bueavov (200(\ ). (b) hadronic emission 



as calculated by Berezh ko &: Volk (20061 ), and (c) IC e mission as cal culated by 



Porter. Moskalenko fc Strong (20061 ). Reproduced from iHinton (20081 ) . 




Figure 8: Multiwavelength views of IC443 (left) and W 28 (right). Molecular 
tracer ^^CO (J=2— t-I ) is shown (cyan contour s) in comparis on to TeV data (colo r 
scale) from H.E.S.S. (|Aharonian et al. 2008el ) and MAGIC (| Albert et al. 2007al ). 
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Figure 9: Five 7-ray PWN candidates in X-rays (left) and TeV 7-rays (right), a) 
Vela X, b) MSH 15-52, c) the K3 and Rabbit PWNe in the Kookaburra Nebula, 
and d) G18. 0-0.7 / HESS J1825-137. Th e 7-ray images are ah made using 



H.E.S.S., see Aharonian et al. (2005a . 2006d 31a) . Publicly available X-ray data 



have been reprocessed to produce the X-ray images: a) ROSAT, b) Chandra, c) 
XMM and Chandra (white inset), d) XMM. The positions of the associated radio 
pulsars are shown with crosses. The white scale bars are 0.5° long. 
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Figure 10: SED for HESS J1640— 465, showing hypothetical contributions from 
young/old electrons. Data for HESS J1825— 1 37 are shown for comparison (gray 
regions). Reproduced from Funk et al. (200?! ). 
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Figure 12: Phase-folded light curv e and spectral index var iations for the binary 
system LS 5039. Reproduced from Aharonian et al. (2006al ). 
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Figure 13: Selected unidentified 7-ray sources as seen by H.E.S.S. The Galactic 
plane is shown as a dashed line. The smoothed PSF is ind icated by a circle at 
the bottom-left of each image. See Aharonian et al. for details. 
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Figure 14: VHE light curve of PKS 2155-304 during the July 2006 flare, 
i n two energy bands: 200-800 GeV (bottom) and above 800 GeV (top) 
( Aharonian et al. 2008al ). The light curve is sampled in two- minute intervals 
around each point. 
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Figure 15: Spectral energy distribution of Mrk 501 in two different states 
( Katarzviiski. Sol fc Kus 200 ll ). The two dominant peaks are interpreted as 
synchrotron and synchrotron self Compton (SSC) emission of electrons (dashed 
lines). Modeling of the SED at lower frequencies adds contributions where ex- 
ternal photons are Compton-scattered (dotted lines) and emission by the host 
galaxy (thin full line). 




Figure 16: EBL limits obtained from VHE 7-ray spectra using plausible assump- 
tions for intrinsic spectra (solid black line), compared to lower limits from direct 
observations — potentially hampered by incomplete subtraction of foregr ound 
emission — and including recent EBL models. From iMazin Raue 12003); see 
there for details and references. 
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Figure 17: Spectral energy distribution from DM annihilation compared to the 
spectrum of 7-rays from the Galactic center as measured using H.E.S.S. Dashed 
line: Annihilation of "typical" 14 TeV WIMPs, solid line: mix of 70% bb and 
30% r+r~ modes of 10 T eV WIMPs, dotted: annihilation of 5 TeV KK particles 
(jAharonian et al. 2006hl ). 
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Figure 18: 7-ray s pectrum resulting from DM annihilati on in one of the bench- 
mark scenarios of ( Bringmann. Bergstrom &: Edsio 20081 ). with strong contribu- 
tions from internal bremsstrahlung. 



